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PRESENTACION 


Este documento de Evaluacion del Potencial Eolico y Solar en El Salvador, se enmarca dentro de los esfuerzos 
que esta desarrollando el Ministerio de Medio Ambiente y Recursos Naturales (MARN), con el apoyo directo de 
la Universidad Centroamericana “Jose Simeon Canas” (UCA), para el cumplimiento de los compromisos 
adoptados por el pais dentro de la Convencion Marco de las Naciones Unidas sobre el Cambio Climatico, 
encaminados a la promocion de medidas, programas y proyectos que reduzcan emisiones de gases de efecto 
invernadero. 

Esta publicacion, es uno de los resultados concretos que se han logrado a traves del proyecto regional 
“Evaluacion del Potencial de Energfa Eolica y Solar” (siglas en ingles SWERA), el cual ha sido coordinado por el 
Programa de las Naciones Unidas para el Medio Ambiente (PNUMA) y apoyado financieramente por el Fondo 
para el Medio Ambiente Mundial (GEF). 

El proposito de formular, editar y publicar este documento es contribuir a superar la barrera de ausencia y 
accesibilidad de informacion confiable, sobre los recursos de energfa solar y eolica en El Salvador; con lo cual 
se espera apoyar la consecucion de los siguientes objetivos especfficos: 

i) Reducir la incertidumbre asociada a las decisiones de inversion y desarrollo para proyectos de 
energfa solar y de viento. 

ii) Evidenciar la abundancia de radiacion solar con la que cuenta la mayor parte del territorio nacional, 
durante los diferentes meses del ano, y potenciar su aprovechamiento. 

iii) Facilitar el diseno tecnico de proyectos de energfa solar. 

iv) Facilitar la identificacion de las zonas geograficas del pais que cuentan con el mayor potencial 
eolico, sirviendo de base para que otras entidades y empresas interesadas evaluen la factibilidad 
de ejecutar en zonas especfficas, tanto proyectos a pequena escala como proyectos comerciales 
de mayor escala. 

v) Contribuir a la reduccion de emisiones de gases de efecto invernadero, a traves de la promocion de 
fuentes renovables de energfa que eviten la quema de combustibles fosiles. 

vi) Fortalecer la estrategia de promocion de las energfas renovables ejecutada por el MARN en 
coordination con el Ministerio de Economfa. 

El equipo de trabajo de la UCA que participo directamente en la elaboration de este documento, estuvo 
conformado por investigadores del Departamento de Ciencias Energeticas y del Laboratorio de Sistemas de 
Informacion Geografica de la Maestrfa en Gestion del Medio Ambiente (MAGMA/UCA). 
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1. Introduccion 


En este documento se presentan las estimaciones de mapas de radiacion solar, brillo solar (horas de luz 
solar) y de potencial eolico para El Salvador. Los mapas de tipo solar han sido generados utilizando 
tecnologia de sistemas de informacion geografica y, en el caso del viento se ha utilizado el software 
WAsP. Para el caso de potencial solar (brillo solar y radiacion solar) se ha utilizado la tecnica de 
interpolation de Kriging. Sus estimaciones se han realizado a partir de mediciones en campo en 22 
estaciones meteorologicas en las cuales se utilizo el modelo de Angstrom y datos de clima (zonas de vida 
de Holdrige) para estimar datos de radiacion solar en estaciones con mediciones unicamente de brillo 
solar. En el caso de potencial eolico las velocidades promedio, direcciones predominantes del viento, los 
calculos de densidades de potencia y estudios de modelaje eolico se han realizado con registros de 18 
estaciones meteorologicas. Durante el desarrollo del proyecto SWERA se han actualizado los datos 
digitales de potencial eolico y solar del pais disponibles en la division de meteorologia del MARN/SNET, 
los avances alcanzados en este proceso se presentan en la section 2 de este documento. 


1.1 Ubicacion geografica v clima de El Salvador. 


La Republica de El Salvador esta localizada en la parte noroeste de Centroamerica, limita al norte y este 
con Honduras, en el extremo sureste con el golfo de Fonseca, al sur con el oceano pacifico, y al oeste y 
noroeste con Guatemala. Geograficamente esta ubicado entre las siguientes coordenadas: 

Latitud 13 s 9'N / 14 2 26'N 
Longitud 87 2 43’0 / 90 s O8’0 

El Salvador es el pais mas pequeno y mas densamente poblado de Centroamerica, tiene su meridiano 
estandar en 89 2 0 (ver Figura 1-1 (a)). Su superficie es de 21,041 Km 2 . El pais esta formado por una 
meseta central recortada por valles fluviales, entre los que destaca el valle del rio Lempa, su cause fluvial 
de mayor longitud. Esta meseta esta flanqueada por dos cadenas montanosas de origen volcanico: la 
primera corre paralela al oceano pacifico y la segunda constituye la frontera con Honduras. Posee una 
estrecha llanura costera en el pacifico de menos de 24 Km de anchura. En el extremo norte se localizan 
altas y escarpadas montanas. La meseta central tiene una altitud promedio de 600 m. El punto mas alto 
del pais es el volcan Santa Ana, con 2,385 m, situado en la cadena costera (ver Figura 1-1 (b)). 

Desde el punto de vista climatologico, El Salvador se encuentra situado en la parte exterior del Cinturon 
Climatico de los Tropicos, caracterizado por tener unas propiedades termicas casi constantes durante 
todo el ano y dos estaciones muy marcadas, una lluviosa y otra seca. Las precipitaciones se concentran 
casi exclusivamente durante la estacion lluviosa, si bien, dentro de ella pueden aparecer grandes 
oscilaciones. Sus valores maximos suelen darse unas semanas despues del paso del sol por su cenit. 

La estacion seca se presenta durante el semestre invernal que comprende los meses de noviembre a 
abril, presentandose las maximas temperaturas al final de la misma, en el mes de abril, antes del 
comienzo de la estacion lluviosa. Otra caracteristica del clima en este pais son los Alisios, vientos 
predominantes de rumbo aproximando NE y los Nodes, que transportan aire frio procedente del Artico 
hacia los tropicos. Se puede observar, tambien, un buen desarrollo del sistema que Neva aire del mar 
hacia tierra y sus zonas inmediatas. El viento se genera por el calentamiento desigual de la superficie de 
la tierra realizado por el Sol (Tabla 1-2). 
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Relieve de El Salva 


Figura 1-1 (a) Ubicacion geografica de El Salvador; (b) Relieve del pais. 


Con los datos obtenidos a lo largo de mas de 50 anos se han podido establecer unas fechas promedio 
para el comienzo y el final de las dos estaciones, asi como para los intervalos de transicion que hay entre 
ellas (Tabla 1-1). 


Tabla 1-1 Duracion aproximada de las estaciones climaticas en El Salvador. 


Epoca del ano 

Principio 

Final 

Dias 

Estacion seca 

14 noviembre 

19 abril 

157 

Transicion seca-lluviosa 

20 abril 

20 mayo 

31 

Estacion lluviosa 

21 mayo 

16 octubre 

149 

Transicion lluviosa-seca 

17 octubre 

13 noviembre 

28 


Fuente: Caracterizacion de los riesgos geologicos y dimensionamiento de los recursos hidrogeologicos de 
Nejapa realizado por Geologos del Mundo y la UCA en agosto de 2003. 
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Tabla 1-2 Regimen de vientos en El Salvador. 


Tipo de viento 

Descripcion 

Vientos alisios de 
caracter regional 

El Salvador esta sometido a la influencia de los vientos alisios, que constituyen los 
vientos mas constantes de todo el planeta . Estos soplan en la direction NE desde 
las regiones de altas presiones de las latitudes medias hacia la zona de calmas 
ecuatoriales o frente intertropical. Estos vientos afectan con mayor intensidad las 
zonas de mayores alturas del territorio. 

Vientos de 
caracter local: 
brisas mar-tierra 

Poseen dos direcciones, la primera se produce en sentido mar-tierra durante el dia, 
debido al gran calentamiento que se opera en la superficie terrestre, que provoca una 
zona de menor presion. Mientras, sobre la superficie del oceano Pacifico, que se 
encuentra mas fria, se genera una zona de mayor presion, lo que provoca que el aire 
circule en sentido oceano-tierra. Estas brisas marinas pueden llegar a penetrar hasta 
unos 100 Km tierra adentro, llegando hasta el valle medio del Rio Lempa. Por la 
noche se produce el fenomeno contrario, debido a la rapida perdida del calor del 
continente, y el mayor calentamiento de la masa de agua del oceano. 

Vientos de 
caracter local: 
brisas montana- 
valles 

Estos vientos poseen dos direcciones en funcion del momento del dia. Durante el dia 
las zonas mas altas del pais se calientan mas rapidamente que los valles y mesetas, 
por lo que se generan vientos locales que soplan desde los valles y mesetas hacia 
las zonas montahosas (brisas del valle), suavizando las temperaturas diurnas de las 
areas de mayor altura del pais. Cuando llega la noche las zonas altas del pais se 
enfrian rapidamente, mientras que las zonas bajas (valles y mesetas) pierden 
lentamente el calor acumulado durante el dia, por lo que las brisas soplan desde las 
zonas montahosas hacia los fondos de los valles y mesetas (brisas de montana), 
suavizando las temperaturas nocturnas de las zonas bajas. 


Nota: Adaptado de la seccion 1 .5 Diagnostico-Clima del “Plan Nacional de Ordenamiento y Desarrollo 
Territorial”. 


1.2 Antecedentes de mediciones solares v provectos 

solares. 


El principal antecedente de determination del potencial solar en el pais lo constituye el proyecto de 
Solarimetrfa ejecutado por el departamento de Fisica, de la Facultad de Ingenieria y Arquitectura de la 
Universidad de El Salvador (UES). Este proyecto, financiado por la Organizacion de Estados Americanos 
(OEA), fue ejecutado de septiembre de 1986 a enero de 1987 y sus resultados fueron publicados en 
septiembre de 1989. 

La UCA con esta investigacion, enmarcada dentro del proyecto SWERA, actualiza con tecnologia la 
recuperacion de datos digitales, la estimacion del potencial solar del pais y amplia la informacion 
disponible en la parte eolica. Tambien la UCA ha realizado en los ultimos cuatro anos instalaciones 
fotovoltaicas en el area rural (4 en total) y en el campus (principalmente con fines de investigacion). 
Destacan en este esfuerzo el proyecto fotovoltaico de El Alto (financiado por Solidaridad Internacional), 
que consistio en la electrificacion de 27 casas, 1 escuela y 1 iglesia a traves de energia solar y el 
proyecto edificio Martin-Baro ((financiado por UCA) que cubre el 30% de la carga de iluminacion del 
edificio (ver Figura 1-2). El sitio web de la UCA http ://dei . uca.edu .sv/dei/htm l/provectos/au lasc/au lac- 
ano.html?ano=3 presenta datos diarios de mediciones solares en Ifnea, registrados en el campus UCA, 
para los anos 2002 a 2004. 

Los datos recopilados en este estudio provienen del Servicio Meteorologico Nacional, ahora conocido 
como Division de Meteorologia del Servicio Nacional de Estudios Territoriales (SNET). Tambien se 
presentan datos de estaciones meteorologicas de la Fundacion Salvadorena para Investigaciones del 
Cafe (ProCafe) que mide radiacion solar en las principales zonas cafetaleras del pais. Los datos de 
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ProCafe fueron extrafdos de los boletines mensuales publicados por esta fundacion a traves de su 
sistema de informacion meteorologica. Es importante mencionar que debido a los pocos anos de 
mediciones realizados por ProCafe sus datos no fueron utilizados para las estimaciones realizadas de 
potencial solar. 


1 .3 Antecedentes de mediciones de viento en el pafs. 


Existen pocos estudios sobre la determinacion del potencial eolico en el pafs. Uno de ellos lo constituye 
“El viento en Centroamerica” realizado por NRECA en la decada de los anos 80. Las principales 
conclusiones de este estudio para el caso de El Salvador fueron: 

■ Se seleccionaron 3 sitios para realizar mediciones de potencial eolico: zona sur del departamento 
de Santa Ana, al este del lago de Coatepeque; zona sur-oeste del departamento de la Libertad; y 
zona centro-oeste del departamento de San Miguel, limftrofe con el departamento de Usulutan. A 
la fecha de publication del documento se tenfan menos de 6 meses de mediciones con sensores 
de velocidad y de direction a 15 y 30 metros de altura en torres de tipo NRG. 

■ Se obtuvo la informacion de 71 estaciones meteorologicas de la region centroamericana que 
median velocidad del viento. Las alturas se homogenizaron a 30 metros de altura utilizando la 
ecuacion: v, = v 0 * (z/z 0 ) a , en donde a es el coeficiente de rugosidad igual a 1/7. Tambien en 
algunas estaciones se corrigieron los datos por un factor de deterioro del equipo de medicion. 

■ Se encontro que en el pafs existe un patron de viento estacional fuertemente marcado. El primero 
es el patron general de vientos de la region (alisios) que tiene su mayor velocidad entre los 
meses de octubre a febrero. El otro sistema es el de las brisas del mar que se contrapone al 
sistema general y tiene sus mayores velocidades en los meses de abril a septiembre. 

■ El Salvador no cuenta con zonas de altas velocidades de viento ; las mas altas estarfan en la 
parte central hacia el norte 5-6 m/s. El resto del pafs tiene velocidades entre 4-5 m/s con un 
patron constante a lo largo del ano. Aunque no se hace explfcito en el documento, se infiere que 
estas velocidades se obtendrfan a 30 m sobre el nivel de la superficie de la tierra. Tambien 
concluye que con la informacion recabada con las mediciones de NRECA mas la informacion 
meteorologica de las 20 estaciones meteorologicas del servicio meteorologico del pafs, la 
posibilidad de que en El Salvador se pueda instalar un parque eolico es reducida, no asf 
proyectos de bombeo de agua, telecomunicaciones y/o generation aislada donde si muestra un 
gran potencial. 

En la section DIAGNOSTICO-clima del “Plan Nacional de Ordenamiento y Desarrollo Territorial” se 
hace referencia a dos estudios que tratan sobre la determinacion del potencial eolico en El Salvador. Las 
principales conclusiones de cada uno de ellos se presentan a continuation: 

(a) “Extraccion de agua subterranea con energia eolica en El Salvador” realizado en la Universidad 
Politecnica de El Salvador en 1987. 

■ El objetivo de este trabajo fue el de determinar el mapa eolico de El Salvador para aprovechar la 
energfa en la extraccion de aguas subterraneas. 

■ Utilizando la informacion de 20 estaciones meteorologicas con mediciones de velocidades del 
viento determino que las velocidades mas altas de 15 minutos se observan de noviembre a 
marzo con valores de 13.9 m/s en la costa, pudiendose alcanzar los 47.2 m/s en las cimas de las 
montanas, por encima de la cota de los 1 600 msnm. 1 


1 En este proyecto las velocidades instantaneas maximas (rafagas) registradas en las estaciones del pafs estan en el 
rango de 4.5 a 7.8 m/s. 
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Figura 1-2 (a) Instalaciones fotovoltaicas proyecto El Alto (Chalatenango); (b) Instalaciones fotovoltaicas 

edificio Martin-Baro (UCA). 


■ Del analisis de los resultados estadfsticos se infiere que en la mayorfa de las estaciones el 
promedio de la velocidad del viento oscila entre los 2.22 - 3.61 m/s. Las estaciones que 
presentan las mayores velocidades medias anuales se presentan en la tabla 1-3. Las menores 
velocidades se presentaron en las estaciones: Estacion Matriz (1.56 m/s), Santa Cruz Porrillo (1.5 
m/s), Beneficio El Papalon (1.5 m/s), San Andres (1.5 m/s), Nueva Concepcion (1.42 m/s) y 
Cerron Grande (1.39 m/s). 
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Tabla 1-3 Areas que presentan mayor potencial eolico de acuerdo a estudio UPES. 


Codigo 

Nombre y depto. 

Velocidad media 
anual (m/s) 

Densidad de potencia 
media anual (W/m 2 ) 

Densidad de energia 
media anual (KWh/m 2 ) 

A-35 

Cerro Verde (Santa 
Ana) 

4.39 

116 

1040 

A-31 

Planes de 
Montecristo (Santa 
Ana) 

3.14 

70 

629 

Z-4 

La Galera (Morazan) 

2.86 

57 

519 


Fuente: Tornado del “Plan Nacional de Ordenamiento y Desarrollo Territorial”. Velocidades a 10 metros 
sobre el nivel del suelo. 


(b) “Primer plan nacional de desarrollo energetico integrado 1988-2000”, publicado en enero de 
1988 por la Comision Ejecutiva Hidroelectrica del Rio Lempa (CEL). 

■ Las conclusiones son similares a las obtenidos por el estudio de la UPES. La mayorfa de 
estaciones de medicion del viento tienen densidades de energia anual comprendidas entre 100 - 
500 KWh/m 2 y solo tres de las 20 estaciones estan por encima de 500 KWh/m 2 , destacandose 
Cerro Verde con 1040 KWh/m 2 . Resultados similares a los presentados en la tabla 1-3. 

En noviembre de 1994 Alexandre de Lemos Pereira del Centro Brasileiro de Energia Eolica realizo con 
apoyo de la UCA, MARN y la CEL el estudio “Site selection for wind power project in El Salvador”. 
Este proyecto se ejecuto con fondos del proyecto SWERA-UCA y con el apoyo del laboratorio de SIG de 
la UCA. El objetivo de este proyecto fue el de seleccionar varios sitios del pais para realizar mediciones 
de viento, hacer recomendaciones para las mediciones a realizar en cada uno de ellos y definir las 
especificaciones tecnicas para el sistema de adquisicion de datos de una torre de medicion de velocidad 
y direccion de viento. Las principales conclusiones de este estudio fueron: 


■ La metodologia para ubicacion de los sitios consistio en: (a) analisis de la informacion disponible 
de viento, en este caso el mapa eolico generado por NREL para el pais del ano 2004 como parte 
del proyecto SWERA y la informacion meteorologica historica recopilado por la UCA en formato 
digital; (b) Los criterios para seleccion del sitio fueron los siguientes: ubicados en regiones de 
densidades de potencia de al menos 300 W/m 2 (correspondiente a 7 m/s y a una altura de 50 m) 
de acuerdo al mapa de NREL, analisis topografico para identificar caracteristicas que 
contribuyeran a mejorar los recursos de viento, tipo de vegetacion, exclusion de areas protegidas, 
buen acceso y proximidad a la red electrica; (c) Con apoyo de personal de la UCA y CEL se 
visitaron 7 diferentes regiones del pais durante 4 dias: area a lo largo de la frontera con 
Guatemala llamada La Hachadura, area al sur del lago de Coatepeque entre Armenia y el cerro 
San Marcelino, area al noroeste de la ciudad de La Libertad llamada Las Shilas, area al node de 
la ciudad de Metapan arriba de la frontera con Guatemala, area alrededor de las ciudades de La 
Palma y San Ignacio, area entre el rio Lempa, San lldefonso y Santa Clara y area alrededor de la 
villa de Monteca. 

■ Despues de aplicar los criterios anteriores se ubicaron cuatro sitios: Metapan, se localizo una 
pequena cresta en la loma de Santa Cruz, cerca de la ciudad de Metapan a 600 msnm 
(89.477977 W y 14.345140 N) (better recommended)', San Isidro, localizado en un valle 
relativamente piano cerca de la villa de San Isidro a una altitud de 780 msnm (89.559193 W y 
13.793907 N) (good recommmended)', La Hachadura localizada cerca de la frontera con 
Guatemala a lo largo de la carretera CA2 en la costa, a una altitud de 50 msnm (90.082879 W y 
13.853474 N) (best recommended)', y, Monteca localizado a 900 msnm en una planicie cercana a 
la frontera con Honduras con acceso diffcil desde la ciudad de Santa Rosa de Lima (87.85512 W 
y 13.87799 N) (recommended). 


Proyecto SWERA-EL SALVADOR 


1-6 



# 


HrOnAi^rttyi 

HttiruaNilmia 



A la fecha de elaboration de este documento se tiene conocimiento de un estudio de viento realizado 
hace varios anos en el pais por OLADE y datos de viento usados por el consultor colombiano Humberto 
Rodriguez quien participa en el proyecto “Implementation de medidas para el desarrollo de un mercado 
national sostenible para servitios de energta en el area rural” para el Ministerio de Economia. 


1.4 Obietivos principales de esta investiqacion. 


El Objetivo principal de esta investigation es la determination del potencial solar y eolico del pais en base 
a datos registrados en las estaciones meteorologicas de la division de meteorologia del SNET. Dado que 
a la fecha no existla una base de datos digital actualizada, los esfuerzos estuvieron dedicados a cumplir 
los siguientes objetivos especfficos: 

■ Actualizar la base de datos digital de datos de radiacion y brillo solar (horas de luz solar) y datos 
de viento (velocidad, direction y rafaga). Para esto fue necesario realizar un inventario completo 
de la informacion existente en el SNET. 

■ Disenar un sistema informatico que permita en el futuro la actualization periodica de informacion 
registrada por las diferentes estaciones meteorologicas del pais. Este sistema a su vez permite 
generar diferentes consultas y reportes de la informacion disponible. 

■ Investigar los diferentes sistemas de interpolation utilizados en Geoestadfstica (IDW, Spline, 
Kriging, etc) para generar los mapas de potencial solar. Finalmente, se selecciono el metodo de 
Kriging para generar los modelos de superficies. 

■ Utilizar el metodo de Angstrom y el mapa de zonas de vida de Holdrige para estimar datos de 
radiacion solar en estaciones con mediciones de luz solar unicamente. El principal aporte en este 
sentido fue la creation de un modelo cartografico, propio de un sistema de informacion 
geografica (SIG), para generar los mapas de radiacion solar y luz solar. 

■ Utilizar el software WAsP para hacer estimaciones de potencial eolico en el pais y utilizar el mapa 
eolico desarrollado por NREL para el proyecto SWERA. 

■ Utilizar WAsP para modelar los sitios de mayor potencial eolico en el pais. 

En las siguientes secciones se presentan los resultados obtenidos los cuales se espera sean de mucha 
utilidad para investigadores y usuarios de informacion solar y eolica. 
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2. Metodologia de la investigacion. 


Para determinar el potencial de energla solar y eolica del pais los esfuerzos iniciales del proyecto 
SWERA-EI Salvador se concentraron en identificar y completar una base de datos digital del pais sobre 
los parametros climatologicos de sol y viento registrados por la division de meteorologia del Servicio 
Nacional de Estudios Territoriales (SNET) y de otras instituciones nacionales. Este proceso consistio de 
las siguientes actividades: digitalizacion de datos existentes, inventario de registros climatologicos del 
Servicio Nacional de Estudios Territoriales (SNET) y ProCafe, creacion de interfaz grafica para la captura 
de radiacion solar en ArcView, creacion del Sistema Automatizado de Manejo de Datos Meteorologicos 
(Luz, Viento y Radiacion), uso de SIG y tecnicas de interpolacion (geoestadistica) y uso del software 
WAsP para estimaciones de potencial eolico. La metodologia utilizada en esta investigacion se ilustra en 
la figura 2-1 . 


Metodologia utilizada para la determinacion del 

potencial solar y eolico 


Estudio del clima, regimen de vientos, 
determinacion de estaciones 
meteorologicas con mediciones de sol 
y viento. 


Antecedentes Solar: 

Proyecto Solarimetria (LIES, 1986-1989) 

Antecedentes eolicos: 

Proyectos NREC.A, UPES, CEL, OLADE, 
MINES, SWERA 

Estudios cientificos de sol y viento 



Inventario de la 
informacion solar y 
eolica de la division 
de meteorologia del 
SNET: 

Radiacion solar 

Horas de luz solar 
(brillo) 

Viento (direction, 
velocidad y rafaga) 


Recuperation de datos 
en formato digital: 
Sistema automatizado 
de datos 

meteorologicos en VB y 
ArcView 


Actualization de datos del proyecto 
solarimetria con datos de SWERA 



Creacion de base de datos de viento 
horaria para potencial eolico en WasP 


* 

Mapa eolico NREL y uso de WAsP para 
determinacion de sitios de mejor potencial 
eolico, calculo de huellas meteorologicas de 
las estaciones del pais. 


I Modelo de Angstrom 

l 


Tecnicas de Geoestadistica: 
Interpolacion de Kriging para 
determinacion de mapas de 
radiacion y brillo solar. 


Figura 2-1 Metodologia de investigacion. 


2.1 Datos del proyecto Solarimetria de la UES. 

Una primera fuente de informacion fueron las tablas de radiacion solar y brillo solar (horas de luz solar) 
publicadas en un documento de septiembre de 1989 como resultado del proyecto “Solarimetria” 
ejecutado por la UES de 1986 a 1987 con el financiamiento de la OEA. En este documento se tenlan 
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registrados datos de radiacion solar de 10 estaciones meteorologicas y brillo solar (horas de luz solar) en 
21 estaciones. Usando el modelo de Angstrom les fue posible estimar la radiacion solar en estaciones 
que tenian registros de brillo solar unicamente pero no de radiacion solar. De esta forma se obtuvieron 
datos de radiacion solar en 22 estaciones del pais (10 con mediciones reales). En la tabla 2-1 aparecen 
las estaciones con mediciones de radiacion solar y en la tabla 2-2 se indican los anos con mediciones de 
radiacion para cada uno de los meses. De manera similar se presentan las tablas 2-3 y 2-4 para el caso 
de horas de luz solar o brillo solar. 


Tabla 2-1 Registros de radiacion solar recopilados por el proyecto de Solarimetna de la UES (Promedios 

mensuales en KWh/mVdia). 1 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a31 

4,5 

4,8 

4,8 

4,2 

3,5 

3,5 

4,7 

4,0 

3,3 

3,8 

4,2 

4,2 

g13 

3,9 

4,1 

4,5 

4,3 

3,8 

3,7 

4,1 

4,2 

3,3 

3,6 

3,4 

3,5 

g3 

4,8 

5,3 

5,6 

5,5 

5,0 

5,1 

5,4 

5,4 

5,0 

5,0 

4,8 

4,6 

h8 

5,1 

5,5 

5,8 

5,6 

5,2 

4,9 

5,6 

5,3 

4,6 

4,8 

4,9 

4,8 

n15 

4,8 

5,2 

5,5 

5,3 

4,8 

5,0 

5,3 

5,2 

4,8 

4,7 

4,6 

4,6 

s27 

5,0 

5,5 

6,1 

5,8 

5,2 

5,2 

5,9 

5,6 

4,9 

4,8 

5,0 

4,8 

s5 

4,8 

5,2 

5,4 

5,4 

4,7 

5,0 

5,5 

5,1 

4,6 

4,6 

4,7 

4,7 

ul 1 

4,9 

5,4 

5,5 

5,2 

4,8 

4,8 

5,4 

5,4 

5,0 

4,8 

4,8 

4,7 

v13 

4,6 

5,2 

5,3 

5,1 

4,8 

4,7 

5,1 

5,1 

4,6 

4,5 

4,6 

4,4 

z4 

4,3 

4,8 

5,1 

4,5 

4,1 

4,0 

4,5 

4,3 

3,5 

3,7 

3,7 

3,9 


Tabla 2.2 Numero de anos con registros de radiacion solar para cada una de las estaciones (Proyecto 

Solarimetna). 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a31 

5 

5 

5 

6 

6 

6 

6 

6 

5 

6 

7 

7 

g13 

6 

6 

6 

6 

5 

5 

4 

4 

4 

5 

6 

6 

g3 

14 

14 

14 

14 

14 

14 

13 

13 

12 

12 

13 

13 

h8 

13 

13 

13 

13 

14 

14 

13 

13 

13 

13 

13 

13 

n15 

13 

13 

13 

13 

12 

12 

14 

13 

13 

13 

14 

13 

s27 

11 

10 

11 

11 

10 

10 

10 

10 

10 

10 

9 

10 

s5 

12 

12 

12 

12 

11 

13 

11 

11 

11 

12 

12 

12 

ul 1 

6 

6 

6 

6 

6 

6 

6 

6 

6 

5 

5 

5 

v13 

8 

8 

8 

8 

8 

8 

8 

9 

9 

8 

8 

8 

z4 

7 

7 

7 

7 

7 

6 

6 

6 

6 

6 

6 

6 


1 Las mediciones de radiacion solar global fueron medidos con Actinografos tipo Robitzch tambien llamado 
Piranografo Bimetalico tipo Robitzch. Los anos de registros corresponden al periodo 1969 a 1983, con variaciones 
en cada una de las estaciones meteorologicas. Segiin el investigador principal del proyecto de Solarimetna, Sra. 
Ethelvina Morillo de Escobar los registros diarios para cada una de las estaciones estudiadas se encuentran en los 
archivos del Departamento de Ffsica de la UES en formato de papel. 
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Tabla 2-3 Registros de brillo solar recopilados por el proyecto de Solarimetria de la UES (Promedios 

mensuales en horas de luz solar). 2 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a18 

7,0 

7,9 

8,3 

7,5 

5,5 

4,7 

6,7 

6,0 

4,5 

4,6 

5,5 

5,9 

a31 

8,8 

8,5 

8,1 

6,7 

5,3 

4,6 

7,0 

5,9 

4,2 

6,0 

7,7 

8,3 

a35 

7,6 

7,3 

7,0 

5,7 

4,9 

4,6 

7,7 

7,0 

4,5 

5,1 

7,0 

7,4 

blO 

8,9 

8,9 

9,5 

8,3 

7,1 

6,7 

8,1 

7,9 

6,4 

7,3 

8,3 

8,8 

P3 

9,1 

9,6 

9,2 

8,5 

7,5 

6,8 

7,9 

7,4 

6,3 

7,3 

8,2 

8,1 

h8 

9,5 

9,7 

9,6 

9,0 

7,5 

6,7 

8,3 

7,4 

5,9 

7,1 

8,6 

9,4 

14 

9,1 

9,1 

8,7 

7,9 

6,9 

6,3 

8,1 

7,4 

5,8 

6,6 

8,2 

8,9 

18 

9,8 

9,8 

9,5 

8,1 

6,7 

6,0 

8,0 

7,4 

5,9 

6,6 

8,6 

9,5 

m18 

9,6 

10,0 

9,7 

8,6 

7,3 

7,0 

8,3 

7,6 

6,3 

7,0 

8,4 

9,2 

m23 

8,0 

8,6 

8,4 

6,2 

4,5 

4,7 

7,2 

6,6 

5,3 

5,5 

6,9 

7,0 

m6 

9,4 

9,8 

9,7 

8,9 

7,7 

7,3 

8,6 

8,6 

7,2 

7,8 

8,5 

9,3 

n15 

9,8 

10,0 

9,7 

8,9 

7,4 

7,4 

8,8 

8,5 

7,2 

8,1 

9,1 

9,7 

s27 

8,9 

9,6 

9,6 

8,6 

6,8 

6,3 

8,3 

7,9 

6,2 

7,0 

8,5 

8,5 

s5 

9,0 

8,7 

9,1 

7,4 

5,5 

5,5 

7,6 

7,3 

5,6 

7,2 

8,3 

8,9 

t6 

9,8 

9,9 

9,6 

8,6 

7,3 

6,8 

8,2 

7,7 

6,7 

7,3 

8,9 

9,4 

u 1 1 

9,7 

10,0 

9,7 

8,7 

7,0 

6,8 

8,1 

8,2 

7,0 

7,3 

8,9 

9,5 

u6 

9,2 

9,4 

9,1 

7,9 

6,4 

5,6 

7,8 

7,4 

5,5 

6,3 

8,0 

9,0 

v13 

9,8 

10,0 

9,5 

8,5 

6,9 

6,5 

8,1 

8,0 

6,6 

7,2 

8,9 

9,7 

v6 

9,5 

9,8 

9,5 

8,5 

7,4 

6,9 

8,3 

8,2 

6,8 

7,4 

8,5 

9,4 

z2 

9,4 

9,8 

9,4 

8,5 

7,0 

7,0 

8,2 

8,0 

6,9 

7,6 

8,6 

9,2 

z4 

7,0 

7,9 

7,9 

7,1 

5,6 

5,1 

7,1 

6,2 

4,6 

4,7 

5,7 

6,2 


Tabla 2-4 Numero de anos con registros de brillo solar para cada una de las estaciones (Proyecto 

Solarimetria). 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a18 

15 

17 

17 

17 

17 

17 

17 

17 

17 

17 

17 

17 

a31 

5 

5 

5 

6 

6 

6 

6 

6 

5 

6 

7 

7 

a35 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

blO 

9 

9 

9 

9 

9 

9 

9 

9 

7 

8 

10 

10 

93 

15 

15 

15 

15 

14 

14 

14 

14 

14 

14 

14 

14 

h8 

10 

10 

10 

10 

11 

11 

11 

11 

11 

11 

11 

11 

14 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

19 

18 

23 

24 

23 

23 

23 

24 

23 

24 

24 

24 

24 

24 

m18 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

m23 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

m6 

21 

21 

21 

21 

21 

21 

21 

21 

21 

20 

21 

22 

n15 

13 

13 

13 

13 

12 

12 

13 

13 

13 

13 

13 

13 

s27 

11 

10 

11 

11 

10 

10 

10 

10 

10 

10 

9 

10 

s5 

5 

6 

6 

6 

6 

7 

6 

6 

6 

6 

7 

7 

t6 

16 

17 

17 

17 

17 

16 

17 

17 

17 

16 

16 

16 

ul 1 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 

10 


2 Las mediciones de brillo solar global fueron medidos con Heliografos tipo Campbell Stokes. Los anos de registros 
corresponden al periodo 1957 a 1984, con variaciones en cada una de las estaciones meteorologicas. 
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Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

u6 

25 

25 

25 

24 

25 

25 

24 

25 

25 

25 

26 

25 

v13 

8 

8 

8 

8 

8 

8 

8 

9 

9 

8 

8 

8 

v6 

16 

16 

16 

17 

17 

17 

17 

17 

17 

17 

17 

17 

z2 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

14 

z4 

8 

8 

8 

8 

9 

9 

9 

9 

9 

9 

9 

9 


2.2 Inventario de informacion solar y eolica de la division de 

meteorologia del (SNET). 

Uno de los objetivos principales del proyecto SWERA fue la recuperacion de informacion solar y eolica de 
los anos posteriores al estudio de la UES (del 984 a 2004). Para este periodo, afectado por la guerra civil 
del pais (1980 a 1992), existfan registros sin procesar; sin embargo, no se encontraban en formato 
digital. La division de meteorologia se limito a archivar unicamente la informacion. Para la recuperacion 
de la informacion la UCA realizo, en primera instancia, los inventarios de la informacion disponible. Estos 
resultados se presentan en las tablas 2-5 a 2- 8. 

Radiacion Solar. 


Para el caso de la Radiacion Solar esta se captura en bandas que contienen siete dias de lectura en 
periodos de 12 horas por dfa, de 6:00 a.m. a 6:00p.m., en la figura 2-1 se muestra un ejemplo de estas 
bandas y en la tabla 2-5 se muestra el inventario total de bandas (2,663). 



Figura 2-1 (a) Ejemplo de una banda de radiacion solar, (b) Metodo del planimetro para recuperacion de 

la informacion de radiacion para el proyecto de Solarimetrla. 

Tabla 2-5 Inventario de informacion de radiacion solar a recuperar en formato digital proyecto SWERA (en 

numero de bandas). 


Codigo 

Estacion 

Departamento 

periodo 

Sub-total 

A31 

Montecristo 

Santa Ana 

1984-1985, 2000, 2002 

120 

H8 

Ahuachapan 

Ahuachapan 

1984-1990, 1992 

255 

G3 

Nueva Concepcion 

Chalatenango 

1 984-2000, 2002 

882 

G13 

Las Pilas 

Chalatenango 

1984-1987, 1989-2000, 2002 

602 

S5 

Observatorio 

San Salvador 

1984-1986, 1996-1998 

215 

S27 

Estacion Matriz 

San Salvador 

1984 

26 

Z4 

La Galera 

Morazan 

1984 

7 

N15 

La Union 

La Union 

1 984, 1 986-2000, 2002 

556 


Total 

2663 
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Brillo Solar (Horas de luz solar) 


Para el parametro de luz solar existen tres tipos de bandas: la curva larga, corta y recta (ver Figura 2-2). 
Cada una de estas bandas se coloca en el heliografo en distintas epocas de ano seleccionadas de 
acuerdo al desplazamiento y orientacion del sol. Cada una de las bandas esta dividida de manera horaria 
en el rango de las 6 a las 18 horas. Para digitalizar se tienen 84,809 bandas. 


Tabla 2-6 Inventario de informacion de luz solar a recuperar en formato digital. 


Codigo 

Estacion 

Departamento 

Anos 

Subtotal 

A18 

Los Andes 

Santa Ana 

1984, 1985, 1986, 1987, 1988, 
1989, 1990, 1991, 1992, 1994, 
1995, 1996, 1997, 1998, 1999, 
2000, 2001, 2002 

4915 

A31 

Los Planes de Montecristo 

Santa Ana 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

6398 

A35 

Cerro Verde 

Santa Ana 

1983, 1984, 1985, 1986, 1987, 
1988 

1554 

A12 

Santa Ana El Palmar 

Santa Ana 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 2002 

4881 


Santa Ana 

Santa Ana 

1986 

29 

HI 

Ahuachapan 

Ahuachapan 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

6623 

H14 

La Hachadura 

Ahuachapan 

1996, 1997, 1998, 1999, 2000, 
2001,2002 

2341 

H99 

San Benito El imposible 

Ahuachapan 

1995, 1996 

577 

T6 

Acajutla 

Sonsonate 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995 

2898 

L31 

Santa Tecla 

La Libertad 

1983, 1984, 1985, 1986, 1987, 
1990, 1991, 1992, 1993, 1994, 
1995,1996, 1997, 1998 

3942 

L4 

San Andres 

La Libertad 

1983, 1984, 1985, 1986, 1987, 
1990, 1991, 1992, 1993, 1994, 
1995, 1996, 1997, 1998, 1999, 
2000, 2001,2002 

5603 

S5 

Observatorio 

San Salvador 

1983, 1984, 1985, 1986, 1987, 
1988, 1989 

1363 

S27 

Estacion matriz(EI Matazano) 

San Salvador 

1983, 1984, 1995,1996, 1997, 
1998 

2600 

G3 

Nueva Concepcion 

Chalatenango 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

5724 
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Codigo 

Estacion 

Departamento 

Anos 

Subtotal 

V6 

Santa Cruz Porrillo 

San Vicente 

1983, 1984, 1985, 1986, 1987, 
1990, 1991, 1992, 1993, 1994, 
1995, 1996 

3772 

BIO 

Cerron Grande 

Cabanas 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

6122 

U11 

Beneficio La Carrera 

Usulutan 

1983, 1984, 1985, 1986 

761 

U6 

Santiago de Maria 

Usulutan 

1983, 1984, 1985, 1986, 1987, 
1990, 1991, 1992, 1993, 1994, 
1995, 1996, 1997, 1998 

4739 

U99 

Laguna de Alegria 

Usulutan 

1998, 1999, 2000, 2001,2002 

866 

M6 

El Papalon 

San Miguel 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

6592 

Z2 

San Francisco Gotera 

Morazan 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

6432 

Z4 

La Galera 

Morazan 

1983, 1984 

231 

M23 

Cerro Cacahuatique 

San Miguel 

1983 

23 

N3 

La Union 

La Union 

1983, 1984, 1985, 1986, 1987, 
1988, 1989, 1990, 1991, 1992, 
1993, 1994, 1995, 1996, 1997, 
1998, 1999, 2000, 2001,2002 

5823 




Total 

84,809 



Figura 2-2 Tres tipos de bandas para medir brillo solar (horas diarias de luz solar). 
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Viento 


Para el parametro viento, los registros se realizan en dos formatos: (a) bandas que son reemplazadas 
diariamente cubriendo periodos de 24 horas, de las cuales se puede extraer, al post-procesarla, datos 
horarios de: rumbo, velocidad y rafaga 3 (23,001 bandas diarias) y (b) Cuadros de registros mensuales 
con informacion de viento (2,558 cuadros). La figura 2-3 muestra las bandas diarias y los cuadros 
mensuales. Los inventarios de viento en bandas y cuadros mensuales se muestran en las tablas 2-7 y 2- 
8 . 
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Figura 2-3 (a) Bandas para registro diario de viento; (b) Cuadro mensual de mediciones de viento (se 
registra velocidad, direccion y rafaga en forma horaria para cada dfa del mes). 

Tabla 2-7 Inventario total de bandas diarias de viento. 


Al 2 

Santa Ana el Palmar 

Santa Ana 

1984-1986, 1994-2000 

2767 

A31 

H8 

T6 

G3 

Planes de Montecristo 

Ahuachapan 

Acajutla 

Nueva concepcion 

Santa Ana 

Ahuachapan 

Sonsonate 

Chalatenango 

1983- 1985, 2000-2002 

1984- 1986, 1994-1996, 1998 
2001 

1986, 1991-2002 
1984-1986 

1472 

2377 

4194 

809 

L4 

San Andres 

La Libertad 

1984, 1990-19997, 1999- 
2002 

4361 

S10 

V6 

llopango 

Santa Cruz Porrillo 

San Salvador 
San Vicente 

1 984-2002 
r 1984 

6221 

N3 

La Union 

La Union 

1984-1986 

707 




Total 

23001 







Rafaga es la velocidad instantanea del viento de mayor magnitud en un dia determinado. 
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Tabla 2-8 Inventario de cuadros mensuales de registros de viento. 


Codigo 

Estacion 

Departamento 

Periodo 

No Cuadros 

A-12 

Sta. Ana El Palmar 

Santa Ana 

1959-1986 

320 

A-31 

Planes de Montecristo 

Santa Ana 

1976-1981 

72 

A-35 

Cerro Verde 

Santa Ana 

1975-1981, 1983 

94 

H-8 

Ahuachapan 

Ahuachapan 

1969-1984 

184 

H-13 

Apaneca Sta. Leticia 

Ahuachapan 

1973-1980, 1982-1984 

128 

T-6 

Acajutla 

Sonsonate 

1954-1956, 1962-1964, 1966- 
1968, 1970-1974, 1976 

173 

L-4 

San Andres 

La Libertad 

1958-1962, 1971-1974, 1976- 
1984 

216 

S-10 

Aeropuerto-llopango 

San Salvador 

1965-1969, 1976-1979, 1983- 
1984, 1986 

132 

S-5 

S.S. Observatorio 

San Salvador 

1955-1974, 1976-1978, 1986 

287 

S-4 

S.S. ITIC 

San Salvador 

1952-1956 

38 

L-18 

El Boqueron 

San Salvador 

1973-1974, 1976-1985 

140 

S-27 

Estacion Matriz 

San Salvador 

1973-1974, 1976-1984 

126 

V-6 

Sta. Cruz Porrillo 

San Vicente 

1969-1982, 1984 

180 

V-13 

Apastepeque 

San Vicente 

1973-1980 

96 

P-11 

Hacienda Astoria 

La Paz 

1970-1975 

120 

B-10 

Cerron Grande 

Cabanas 

1974-1981 

83 

G-3 

Nueva Concepcion 

Chalatenango 

1970-1974, 1976-1985 

169 


Total 

2558 


2.3 Procedimiento para recuperacion de informacion solar y eolica. 


Una vez concluido el inventario de la informacion solar y eolica a recuperar se inicio la creacion del 
sistema automatizado para el manejo de datos meteorologicos en Visual Basic, que servira en el futuro 
como medio de almacenaje y consulta de informacion una vez concluido el proyecto. Dada la cantidad de 
informacion a procesar y el tiempo que se requeria para hacerlo en la forma tradicional, se opto por 
desarrollar el sistema informatico con una interfaz grafica en ArcView que permitiera la captura de bandas 
escaneadas de radiacion solar. A continuacion se detalla el procedimiento para digitalizar cada una de las 
variables meteorologicas de interes a este proyecto. 

Radiacion Solar. 


La interfaz grafica tenia la funcion de calcular los valores meteorologicos de radiacion a partir de las 
bandas analogicas escaneadas. Fue creada con herramientas del SIG ArcView y fue considerada de 
caracter temporal mientras esta tarea se incorporaba al sistema en Visual Basic. 

La interfaz desarrollada permite manejar los archivos escaneados, calibrar las bandas, cambiar las 
unidades de pixeles de una imagen escaneada a unidades reales, dibujar puntos, arcos y poligonos o 
cortar poligonos. Los poligonos posteriormente se convierten a radiacion solar. Los datos calculados se 
almacenaron provisionalmente en Access. A continuacion se describe el proceso para “leer” las bandas 
de radiacion (ver figura 2-4 y 2-5): 

1) Escaneo de las bandas . este paso tiene un doble objetivo, el primero es crear un banco digital de 
imagenes(en formato JPG para generar archivos de calidad con poco tamano) de las 
correspondientes bandas de registro para tener un respaldo si las bandas en papel se deterioran 
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o extravfan; y el segundo, consiste en utilizar esta imagen para la digitalizacion en pantalla de las 
areas a partir de las cuales se calcula la radiacion. Cada imagen se identifica a traves del codigo 
de la estacion mas las fechas de inicio y fin de la lectura. 

2) Digitalizacion en pantalla , una vez se tiene la imagen, el tecnico la carga dentro del Arcview para 
proceder a dibujar las areas que no son mas que la cantidad de radiacion dividida en los periodos 
de manana y tarde. Este proceso genera un Archivo de formato SHP, propio del programa, el 
cual calcula de forma automatica el area bajo la curva del registro. Este dato es almacenado 
directamente en una base de datos para su consulta posterior. 


JLicVIcw CS Vciiiin 3.1 



Figura 2-4 Banda de radiacion solar escaneada e incorporada en ArcView para su digitalizacion. 



Figura 2-5 (a) Generation del banco de imagenes de bandas de radiacion; (b) Calculo de la radiacion 
solar en el software ArcView. 
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brillo solar. 


El procedimiento para digitalizar las bandas de luz solar es el siguiente: cada espacio de hora se divide 
mentalmente en 10 zonas para facilitar la estimacion de la cantidad de tiempo que el sol perforo la banda, 
este tiempo se registra de manera horaria en el computador. Luego, cada una de las estimaciones 
horarias son acumuladas por el software de captura de datos para estimar el numero de horas de luz 
solar en ese dia en particular. Antes de leer la banda se crea un registro en la base de datos de la fecha 
del dla y de la estacion (ver Figura 2-6). 



Figura 2-6 Proceso de captura de luz solar en formato digital. 


Viento. 


El procedimiento para digitalizar la informacion de viento sigue dos procesos: lectura de bandas y lectura 
de cuadros. 

Lectura de banda: se escanea la banda en dos partes (debido a su larga longitud) luego se unen 
mediante software y se identifica el archivo imagen generado con la fecha de registro y el codigo 
de la estacion. Luego se captura en ArcView la imagen escaneada y se capturan tres tipos de 
informacion: direccion del viento, velocidad y rafaga. 

Lectura de cuadro: se digitaliza la informacion de los tres parametros (direccion del viento, 
velocidad y rafaga) los cuales ya estan disponibles en horas, esto debido a que el cuadro es el 
resultado mensual de las lecturas diarias de las bandas para cada una de las estaciones (ver 
Figura 2-7). 
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Figura 2-7 Proceso de captura en formato digital de los cuadros mensuales de viento. 


2.3.1 Sistema Automatizado de Manejo de Datos Meteorologicos 

(Brillo solar, radiacion solar y Viento. 


El Sistema de Captura de Datos, es un software desarrollado a la medida para la recuperacion de datos 
historicos de Viento, Radiacion y Brillo Solar, y tiene como objetivo principal brindar una herramienta 
automatizada para el uso y administracion de los datos de potencial solar y eolico del pais. Este software 
ha sido disenado de la manera mas sencilla, con una interfaz orientada al usuario, permitiendo que la 
aplicacion tenga una interaccion agradable, facil de entender y manejar por cualquier tipo de usuario con 
conocimientos basicos de computacion. 

Mediante esta herramienta se obtendran mayores beneficios en cuanto a las necesidades de informacion 
que surjan, como por ejemplo elaboracion de informes o reportes, ofreciendo una informacion veraz y 
confiable, asf como tambien un analisis de la informacion. 

El Sistema desarrollado facilita el manejo de datos y la obtencion de la informacion a la cual se le da un 
tratamiento mas agil y efectivo. Ha sido desarrollado en el lenguaje de Programacion de Visual Basic 6.0, 
utilizando una base de datos en Access 2000. 

Ventajas del sistema 

• Al instalarse una base de datos centralizada y con conexiones de red, se garantiza la unidad de 
la informacion, facilitando las labores de mantenimiento, respaldo y consulta. 

• El Sistema mecanizado y enlazado por medio de una red con estaciones de usuarios, permite la 
facilidad de consulta, ast, si se necesita algun reporte o informacion procedente de una consulta, 
la puede obtener de forma directa sin necesidad de recurrir a recursos adicionales. 

• El sistema desarrollado brinda el apoyo necesario para el manejo productivo de la informacion y 
la obtencion de la misma, aumentando la capacidad de administrar la informacion. 
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• La informacion puede ser accedida y ordenada segun los criterios que se estimen convenientes, 
evitando asi, la duplicidad de informacion. 


De manera general la estructura del sistema informatico es la siguiente: 

Catalogos 

Esta parte involucra todos los mantenimientos de la aplicacion, algunos de estos mantenimientos 
dependen de otros, como por ejemplo Estaciones y Categorfas o Departamentos y Zona Departamental. 
Cada uno de estos mantenimientos deberan ser llenados antes de utilizar la parte de Rubros . 

a) Categorias 

b) Rumbos 

c) Tipo de Medicion 

d) Zona Departamental 

e) Departamentos 

f) Estaciones 

Rubros 

Esta parte del menu es practicamente donde se ingresan los datos de los tres tipos de Medicion, para 
cada uno de los cuales existe una pantalla diferente. Para los rubros de Viento y Luz Solar, la insercion 
de los datos es practicamente igual, para el rubro de Radiacion el dato se obtendra de una interfaz grafica 
en la cual se escanea la banda y se dibuja sobre ella, realizando una calibracion para obtener un factor, 
el cual se utiliza para el calculo del area. Sin embargo la pantalla para Radiacion es igual que para Viento 
y Luz Solar. De igual manera el rubro de Viento tendra una interfaz grafica, en la cual tambien se 
escaneara la banda pertinente. 

a) Captura de Datos Manual 

1 . Viento 

2. Brillo Solar 

b) Captura de Datos Grafica 

1 . Radiacion 

2. Viento 


Consultas 

a) Consulta de Radiacion 

b) Consulta de Usuarios 

c) Consulta de Estaciones por Categorfa 

Reportes 

a) Mediciones por Estacion 

b) Mediciones por Ano 

c) Medicion por Rango de Anos 

d) Generales 

e) Por Status 

f) Estaciones 

1 . Por Categoria 

2. Por Departamento 

Algunas de las pantallas principales del sistema informatico se presentan en la figura 2- 8 . 
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Figura 2-8 Pantallas del sistema informatico desarrollado para capturar y almacenar datos de potencial 
solar y eolico del pais. 


2.4 Resultados obtenidos del proceso de digitalizacion de la 

informacion en los inventarios. 


A la fecha de elaboracion de este documento el avarice de recuperacion de la informacion digitalizada es 
la siguiente: 

Radiacion solar. 


En el caso de la radiacion solar se recupero en formato digital el 100 % de las bandas (2663). Los 

anos y las estaciones recuperadas se muestran en la tabla 2-9. De las 10 estaciones con mediciones de 
radiacion solar unicamente se han actualizado 8, completando los anos de 1969 a 2002 (en el mejor de 
los casos) para algunas de las estaciones. En la tabla 2-10 se muestran los valores en KWh/m 2 /dia 
actualizados con los datos de Solarimetna y SWERA y en la tabla 2-1 1 se muestra el numero de anos 
con mediciones digitalizadas para cada una de las estaciones. 

Tabla 2-9 Total de bandas digitalizadas de radiacion solar (SWERA). 


Codigo 

Estacion 

Departamento 

periodo 

Sub-total 

A31 

Montecristo 

Santa Ana 

1984-1985, 2000, 2002 

120 

H8 

Ahuachapan 

Ahuachapan 

1984-1990, 1992 

255 

G3 

Nueva Concepcion 

Chalatenango 

1 984-2000, 2002 

882 

G13 

Las Pilas 

Chalatenango 

1984-1987, 1989-2000, 2002 

602 

S5 

Observatorio 

San Salvador 

1984-1986, 1996-1998 

215 

S27 

Estacion Matriz 

San Salvador 

1984 

26 

Z4 

La Galera 

Morazan 

1984 

7 

N15 

La Union 

La Union 

1 984, 1 986-2000, 2002 

556 


Total 

2663 


Proyecto SWERA-EL SALVADOR 


2-13 





Tabla 2-10 Radiacion solar en KWh/m 2 /dia incluyendo datos de Solarimetria y SWERA. 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a31 

4.9 

5.2 

5.0 

4.5 

3.6 

4.1 

5.2 

4.7 

3.8 

4.4 

4.6 

4.6 

g13 

4.6 

5.0 

5.6 

5.3 

5.1 

5.0 

5.5 

5.5 

4.8 

4.7 

4.3 

4.3 

03 

4.9 

5.4 

5.7 

5.5 

5.2 

5.2 

5.4 

5.4 

5.0 

4.9 

4.7 

4.6 

h8 

5.0 

5.4 

5.7 

5.5 

5.0 

4.7 

5.2 

5.0 

4.3 

4.5 

4.6 

4.5 

n15 

4.7 

5.1 

5.4 

5.3 

4.8 

5.0 

5.2 

5.2 

4.8 

4.7 

4.5 

4.4 

s27 

5.1 

5.6 

6.2 

5.9 

5.2 

5.3 

5.9 

5.6 

4.9 

4.8 

5.0 

4.8 

s5 

4.9 

5.4 

5.7 

5.4 

4.9 

5.1 

5.5 

5.2 

4.6 

4.8 

4.8 

4.8 

ul 1 

4.9 

5.4 

5.5 

5.2 

4.8 

4.8 

5.4 

5.4 

5.0 

4.8 

4.8 

4.7 

v13 

4.6 

5.2 

5.3 

5.1 

4.8 

4.7 

5.1 

5.1 

4.6 

4.5 

4.6 

4.4 

z4 

4.4 

4.9 

5.1 

4.8 

4.1 

4.0 

4.5 

4.3 

3.5 

3.7 

3.7 

3.9 


Tabla 2-1 1 Numero de anos de radiacion solar digitalizados para las estaciones meteorologicas (SWERA 

+ Solarimetria). 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a31 

7 

7 

6 

7 

7 

9 

9 

9 

8 

9 

10 

10 

gi3 

19 

18 

19 

18 

18 

18 

18 

17 

17 

15 

16 

21 

g3 

30 

31 

30 

30 

30 

29 

28 

30 

30 

28 

28 

30 

h8 

17 

18 

18 

18 

19 

20 

19 

19 

19 

19 

19 

20 

n15 

23 

23 

23 

23 

22 

24 

25 

26 

26 

25 

25 

24 

s27 

12 

11 

12 

12 

11 

11 

10 

10 

10 

10 

9 

10 

s5 

17 

17 

17 

18 

17 

19 

16 

15 

15 

16 

16 

16 

ul 1 

6 

6 

6 

6 

6 

6 

6 

6 

6 

5 

5 

5 

v13 

8 

8 

8 

8 

8 

8 

8 

9 

9 

8 

8 

8 

z4 

8 

8 

7 

8 

7 

6 

6 

6 

6 

6 

6 

6 


Para propositos de consulta se incluyen los datos de radiacion solar registrados en estaciones 
meteorologicas de Procafe, aun cuando estos datos no se incluyeron para las estimaciones de potencial 
solar. 


Tabla 2-10 (a) Datos de radiacion solar KWh/nr/dfa registrados por estaciones de ProCafe. 


Estacion 

Latitud 

Longitud 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

Ciudad Barrios 

13,44 

-88,3 

6,6 

5,7 

5,0 


5,5 

5,7 

5,6 

6,9 

6,2 

5,4 

7,5 

7,8 

Chalchuapa 

13,93 

-89,67 

4,3 

4,8 

4,5 

6,4 

5,3 

5,2 

5,8 

6,1 

4,8 

5,0 

5,6 

4,4 

Chaparrastique 

13,76 

-88,21 

2,4 

4,6 

4,4 

3,6 

2,4 

3,6 

3,8 

3,2 

2,4 

2,3 

2,4 

2,2 

El Boqueron 

13,76 

-89,29 

5,2 

5,8 

5,2 

7,0 

4,4 

5,0 

5,6 

5,8 

4,7 

4,1 

5,1 

4,9 

Jayaque 

13,71 

-89,45 

5,6 

5,9 

6,2 

7,0 

4,1 

5,1 

5,6 

5,3 

4,3 

4,9 

5,4 

5,8 

Juayua 

13,83 

-89,74 

3,4 

3,9 

3,9 

4,6 

3,2 

3,5 

3,8 

3,7 

2,8 

3,1 

3,3 

3,1 
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Tabla 2-10 (b) Numero de anos de radiacion solar digitalizados para las estaciones meteorologicas de 

ProCafe. 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

Ciudad Barrios 

3 

2 

1 

1 

2 

2 

2 

2 

2 

3 

2 


Chalchuapa 

2 

2 

2 

1 

1 

2 

3 

1 

1 

1 

2 

1 

Chaparrastique 

3 

1 

1 

1 

1 

2 

2 

2 

2 

3 

3 

2 

El Boqueron 

3 

2 

2 

1 

1 

3 

3 

2 

2 

3 

3 

2 

Jayaque 

2 

2 

2 

1 

1 

3 

3 

2 

2 

2 

4 

1 

Juayua 

3 

2 

2 

1 

1 

3 

3 

2 

2 

3 

4 

2 


Brillo solar. 


En el caso del brillo solar, dado que la secuencia de datos de solarimetrfa cubrfa el rango de anos de 
1957 a 1984 para 21 estaciones meteorologicas con datos para los diferentes meses desde 5 anos de 
mediciones hasta 25 anos, se considero que era una buena secuencia de datos para la interpolacion a 
realizar (estimacion del mapa de Luz Solar). Debido a esto el esfuerzo en tiempo de los tecnicos de la 
UCA se dedico a obtener un inventario ordenado de las bandas existentes para las diferentes estaciones 
(84,809 bandas). A la fecha unicamente se han digitalizado los ultimos anos de 17 estaciones, es 
decir, 7,854 bandas que representan un 9.26 % del total (ver tabla 2-12). 

Tabla 2-12 Datos de luz solar digitalizados dentro del proyecto SWERA. 


Codigo 

Estacion 

Departamento 

Anos 

Subtotal 

A18 

Los Andes 

Santa Ana 

2001 , 2002 

508 

A31 

Los Planes de Montecristo 

Santa Ana 

2001, 2002 

695 

A35 

Cerro Verde 

Santa Ana 

1987, 1988 

379 

A12 

Santa Ana, El Palmar 

Santa Ana 

1998, 2002 

432 

H14 

La Hachadura 

Ahuachapan 

2001 , 2002 

726 

H99 

San Benito, El imposible 

Ahuachapan 

1995, 1996 

577 

H8 

Ahuachapan 

Ahuachapan 

2002 

214 

L4 

San Andres 

La Libertad 

2001, 2002 

695 

S27 

Estacion matriz (El Matazano) 

San Salvador 

2000, 2001 

700 

S5 

Observatorio 

San Salvador 

1988, 1989 

224 

BIO 

Cerron Grande 

Cabanas 

2001 , 2002 

698 

U6 

Santiago de Marfa 

Usulutan 

1997, 1998 

392 

U11 

Beneficio La Carrera 

Usulutan 

1985, 1986 

357 

U99 

Laguna de Alegrfa 


2001, 2002 

366 

V6 

Santa Cruz Porrillo 

San Vicente 

1994 

209 

VI 6 

Lempa-Acahuapa 

San Vicente 

2001, 2002 

682 




Total 

7854 
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Viento. 


La informacion recuperada a la fecha es de 719 cuadros mensuales, equivalentes a 20,890 
registros diarios con informacion horaria. La informacion de viento en bandas (23,001) no se han 
digitalizado. 


Tabla 2-13 Cuadros de viento digitalizados en el proyecto SWERA. 


Estacion 

Nombre 

Cuadros 

Dias de medicion 

Anos 

a12 

Santa Ana El Palmar 

68 

2014 

1 980-1 986 

a31 

Los Planes de Montecristo 

36 

1044 

1 980-1 983 

a35 

Cerro Verde 

33 

978 

1980, 1981, 1983 

M23 

Cerro Cacahuatique 

25 

678 

1 980-1 982 

M6 

El Papalon 

48 

1450 

1980-1984, 1986 

N15 

La Union 

22 

667 

1980, 1985 

h8 

Ahuachapan 

50 

1429 

1 980-1 984 

hi 3 

Apaneca Santa Leticia 

35 

997 

1980, 1982-1984 

slO 

Aeropuerto llopango 

53 

1491 

1980-1984, 1986 

s27 

Estacion Matriz 

54 

1599 

1 980-1 984 

14 

San Andres 

54 

1616 

1 980-1 984 

118 

Boqueron 

60 

1756 

1 980-1 985 

blO 

Cerron Grande 

27 

703 

1 980-1 984 

93 

Nueva Concepcion 

61 

1818 

1980-1984, 1986 

v6 

Santa Cruz Porrillo 

48 

1350 

1 980-1 984 

v13 

Apastepeque 

6 

177 

1980 

U11 

Beneficio La Carrera 

17 

505 

1980, 1983, 1984 

Z4 

La Galera 

22 

618 

1 980-1 982 


TOTAL 

719 

20,890 
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3. Aplicacion de Tecnologias SIG al modelo de 

Angstrom 


3.1 Conceptos basicos de Radiacion Solar v Luz solar. 


El Sol v la Tierra 

Se entiende como radiacion solar a la radiacion electromagnetica que proviene del Sol 1 , la cual se 
extiende desde los rayos gamma (de longitud de onda 10' 10 cm e inferiores) hasta las ondas de radio 
(de longitudes de onda 10 +5 cm y superiores). La mayoria de la radiacion solar cae entre las 
longitudes de onda de 0.15 a 120 pm (micrometros), y la radiacion de importancia practica a los 
usuarios de la energia solar cae entre 0.15 y 3.0 pm. Las longitudes de onda de la radiacion visible 
estan entre 0.38 y 0.72 pm (ver Figura 1). La cantidad de energia radiada por el sol se estima a una 
razon de 3.8x1 0 22 kW. 
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Figura 3-1 Espectro Electromagnetico. La radiacion de importancia practica a los usuarios de la 

energia solar cae entre 0.15 y 3.0 pm. 

La Tierra esta a 150 millones de km del Sol y tiene un radio de aproximadamente 6360 km. La Tierra 
gira alrededor del Sol en una orbita eliptica a una velocidad media de 30 km/s y al mismo tiempo rota 
sobre su eje a una velocidad de 0.5 km/s. El eje de rotacion de la Tierra esta inclinado 23. 45 2 con 
respecto a su orbita alrededor del Sol. La distribucion de la radiacion solar sobre la superficie de la 
Tierra y su variation a lo largo de la duration del dia es el resultado de la position inclinada de la 
Tierra, su rotacion diaria y su traslacion anual. Es decir, la radiacion solar recibida en la superficie de 
la Tierra no es constante y para su comprension es necesario conocer los angulos entre el Sol y la 
Tierra. 


1 TA 

El sol es una estrella de radio 0.7 millones de Km, tiene una masa de 2x10 kg y su temperatura en la 
superficie es de aproximadamente de 5760 K. 
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Angulos entre el Sol v la Tierra. 



La posicion de una localidad en la superficie de la Tierra se conoce a traves de las coordenadas de 
latitud y longitud: 

• Latitud (<p) se define como la distancia angular de un punto (en la superficie de la Tierra) 

medida desde el ecuador. Puede ser al node o al sur del ecuador y nos indica que tan lejos 
esta la localidad del ecuador. 

• Longitud (L) o meridiano son semi circulos que pasan por los polos de la Tierra. El meridiano 
cero pasa por Greenwich y por acuerdo internacional se le conoce como meridiano principal. 
La longitud se define como el angulo entre el meridiano principal y el meridiano que pasa a 
traves de una localidad. Localidades al este del Greenwich se definen como grados Este y 
localidades al oeste del Greenwich se definen como grados Oeste. 

Tiempos de reloj (WT). Se basa en longitudes y depende del meridiano estandar para cada pais o 
region. La diferencia en tiempo del meridiano estandar con el Greenwich se conoce como GMT. La 
velocidad de rotacion de la tierra es de 4 minutos por meridiano (o longitud) o de 15 s por hora. 

Tiempo solar (TS). Se basa en el dia solar que se define como el intervalo de tiempo entre dos pasos 
sucesivos del Sol a traves del meridiano del observador. Todos los calculos en energia solar se 
basan unicamente en el tiempo solar. El tiempo solar es diferente del tiempo de reloj por dos razones: 

• El tiempo de reloj depende del meridiano estandar del pais o region y, por lo tanto, existe una 
diferencia entre la longitud local y el meridiano estandar delta en minutos: 

A = 4 (L es -Lj Ec. 3-1 

• Las velocidades de rotacion y angular de la tierra no son constantes. Esta correction se 
conoce como la Ecuacion del Tiempo (EOT). Para calcular EOT se utiliza la siguiente 
ecuacion (en horas): 


EOT = 0.123cos(f + 87)-(l/6) sen [2(e + 10)] E 

s = 0.938[W + 30.3(M - 1)] E 

En donde, N y M representan el numero del dfa en el mes y el numero del mes del ano. 


La diferencia entre el tiempo solar y el tiempo de reloj se relaciona con la EOT de la siguiente 
manera: 


• Para localidades al Este del Greenwich: 
ST = WT + EOT - A 


Ec. 3-4 


• Para localidades al Oeste del Greenwich 

ST = WT + EOT + A Ec. 3-5 

La posicion del sol en el cielo se conoce por los angulos horarios y la declination y estos se 
determinan de la siguiente forma: 
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Angulo horario. 



Es el desplazamiento del sol Este u Oeste del meridiano de una localidad debido a la rotacion 
de la tierra sobre su eje a 15 s por hora. Al medio dla solar (el momento en el cual el sol esta 
sobre el meridiano local) el angulo horario es cero. En la rnahana es negativo, mientras que 
en la tarde, es positivo. El angulo horario se calcula por la siguiente ecuacion: 


a> = 15(57-12) 


Ec. 3-6 


En donde ST es el tiempo solar en horas. 

• Declinacion. 

Es la posicion angular del sol al medio dla solar con respecto al piano del ecuador. Tambien 
se define como el angulo formado por la llnea que se extiende desde el centra del sol al 
centra de la tierra y la proyeccion de esta llnea en el piano ecuatorial de la tierra. La maxima 
declinacion se obtiene en el solsticio de verano (junio 22, 23.45 s ), la minima en el solsticio de 
invierno (diciembre 21 , -23.45 s ) y es igual a cero durante los dos equinoccios (marzo 21 y 
septiembre 22). Esto nos indica que el sol esta directamente sobre nuestras cabezas dos 
veces al ano en localidades entre 23.45 s N y 23.45 s S y el sol nunca esta sobre nuestras 
cabezas a latitudes fuera de este rango. 


8 = 23.45Sen{360(284 + N )/ 365 }grados 


Ec. 3-7 


3.2 Modelo de Angstrom. 

Los datos de radiacion registrados por una estacion meteorologica son la mejor fuente de informacion 
para estimar la radiacion solar incidente promedio en una localidad. Cuando estos datos faltan o no 
hay datos de localidades cercanas, con clima similar, al punto de interes, se utilizan ecuaciones 
emplricas para estimar la radiacion a partir de las horas de luz solar o nubosidad (“cloudiness”). Los 
datos de horas de luz solar promedio (tambien llamado brillo solar en este documento) usualmente 
estan disponibles para mayor numero de estaciones meteorologicas, para el caso de El Salvador. 
Uno de los modelos mas conocidos para estimar datos de radiacion solar es la ecuacion de 
regresion de Angstrom que relaciona, para una localidad de interes, el cociente de la radiacion 
diaria promedio mensual y la radiacion en dla despejado (“clear day radiation)” (H/Hc), con la fraccion 
de horas de luz solar posibles (n/N): 



a +b 


n 

N 


Ec. 3-8 


En donde, 

H : Radiacion diaria promedio mensual en una superficie horizontal. 

H c : Radiacion diario promedio en cielo abierto para la localidad y mes de interes. 

a' y b’: Constantes emplricas. 

n Promedio mensual de horas diarias de luz solar. 

N : Promedio mensual de las horas maximas de luz solar diarias (es decir, la duracion del 

dla del dla promedio del mes). 


Una dificultad de esta ecuacion esta en la ambiguedad de los terminos n/N y H c . El primero es un 
error de tipo instrumental (el registro de luz solar diario esta sujeto a la interpretacion de la persona 
que realiza la lectura en la banda), el segundo, resulta de la incerteza en la definicion de dla 
despejado (“clear day)”. Page y otros autores modificaron la ecuacion original de Angstrom por un 
metodo basado en la radiacion extraterrestre sobre una superficie horizontal en lugar de la radiacion 
diaria promedio en cielo despejado, 
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H , n 
= a + b- 

H o N 

En donde Ho es la radiacion fuera de la atmosfera para la misma localidad promediada en el periodo 
de interes y a y b son las constantes modificadas las cuales, obviamente, dependen de la localidad en 
estudio. 

La metodologia empleada para evaluar la energia solar disponible de las estaciones que poseen 
informacion basica completa (radiacion solar + brillo solar) es la siguiente: 



Ec. 3-9 


Paso 1 . 

Calculo de los promedios mensuales de la duracion astronomica del di'a (longitud del dia u 

horas maximas de luz solar diarias). 


La duracion astronomica del dia N d , o longitud del dia, se calcula mediante: 


2 2 

N d = — co s = — cos -1 (- tan 8 tan <p) 
1 5 1 5 


Ec. 3-10 


en donde: 

° : Declinacion, es la posicion angular del sol al mediodia solar con respecto al piano del 

ecuador. 



^0.006918 - 0.399912 cos y + 0.070257 sin y - 0.006758 cos 2 y^ 
v + 0.000907 senly - 0.002697 cos 3 y + 0.00148 sen3y y 


180 

n 


Ec. 3-1 1 


En donde: 

/17 : es el angulo del dia(en radianes) y se representa por: 


7 = 


2n{d n - 1) 


365 


Ec. 3-12 


Donde d n es el numero del dia (1 < d n < 365). 

^ : Latitud, es el angulo con que se determina el lugar de interes (localidad) sobre la 

tierra con respecto al piano del ecuador. 

Ct>s : Angulo horario con que se oculta el sol. 


El promedio mensual de la duracion astronomica del dia N d puede calcularse, con bastante 
aproximacion, usando en la ecuacion 3-10 el valor de la declinacion para el dia recomendado para 
cada mes. 2 En la tabla 3-1 se muestra el calculo del angulo de declinacion para cada dia 
caracteristico del mes y en la tabla 3-2 se presenta la duracion astronomica del dia para cada mes y 
en cada una de las estaciones consideradas. 


" El dia recomendado de cada mes o dia promedio del mes es el dia en el cual su radiacion extraterrestre es la 
mas cercana al promedio del mes. 
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Tabla 3-1 Calculo del angulo de declinacion 


# 


VErdruL^ 

m v- . ib . 

Hn.^nn^dliiiLEk 



d n 

r 

° (grados) 

3 (rad) 

17 

0.275427945 

-20.90354682 

-0.364835717 

47 

0.791855342 

-12.60892623 

-0.220067278 

75 

1 .273854247 

-2.041903131 

-0.035637933 

105 

1.790281644 

9.480836087 

0.165471806 

135 

2.306709041 

18.67365533 

0.325916769 

162 

2.771493699 

23.0378951 

0.402087122 

198 

3.391206575 

21.34544923 

0.372548369 

228 

3.907633973 

13.98911474 

0.244156112 

258 

4.42406137 

3.34279361 1 

0.058342755 

288 

4.940488767 

-8.217977789 

-0.14343077 

318 

5.456916164 

-18.04105059 

-0.314875733 

344 

5.904486575 

-22.84065258 

-0.398644591 


Tabla 3-2 Duracion astronomica del dfa para cada estacion (en horas). 


Codigo 

Latitud 

Longitud 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

H-8 

13.94 

-89.86 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

A-18 

13.88 

-89.64 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

A-31 

14.4 

-89.36 

11.2 

11.6 

11.9 

12.3 

12.7 

12.8 

12.8 

12.5 

12.1 

11.7 

11.4 

11.2 

A-35 

13.83 

-89.62 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

T-6 

13.57 

-89.83 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

L-4 

13.81 

-89.41 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

L-8 

13.69 

-89.29 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

S-5 

13.7 

-89.2 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

S-27 

13.69 

-89.14 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

G-3 

14.12 

-89.29 

11.3 

11.6 

11.9 

12.3 

12.7 

12.8 

12.8 

12.5 

12.1 

11.7 

11.4 

11.2 

G-13 

14.36 

-89.09 

11.3 

11.6 

11.9 

12.3 

12.7 

12.8 

12.8 

12.5 

12.1 

11.7 

11.4 

11.2 

V-6 

13.44 

-88.8 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

V-13 

13.7 

-88.75 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

U-6 

13.48 

-88.47 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

U-11 

13.33 

-88.52 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

M-6 

13.44 

-88.12 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

M-18 

13.71 

-88.36 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

M-23 

13.76 

-88.23 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

7-2 

13.7 

-88.11 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

Z-4 

14.04 

-88.09 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

N-15 

13.33 

-87.88 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 

B-10 

13.94 

-88.91 

11.3 

11.6 

11.9 

12.3 

12.6 

12.8 

12.7 

12.5 

12.1 

11.7 

11.4 

11.2 
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Paso 2. 


Calculo del promedio mensual de la radiacion solar 

extraterrestre en superficie horizontal. 

La radiacion solar extraterrestre en la superficie horizontal se determina por medio de la siguiente 
ecuacion: 


24(3600) 


0 


l + 0.033cos 


f 360 d n \ 


n 


v 


365 


J 


cosmos S«n», Aen<pse„S 


J Ec. 3-13 


En donde, I sc = Constante solar y su valor es 1353 watt/m 2 y los otros valores fueron definidos 
anteriormente. 

La radiacion extraterrestre, calculada a partir de la ecuacion 3-13, para cada una de las estaciones se 
presenta en la tabla 3-3, de estas 22 estaciones, solo 9 tienen registros de brillo y radiacion como se 
menciono en la seccion 2 de este documento. 


Tabla 3-3 Radiacion extraterrestre para cada estacion (KWh/m 2 /dia) 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a18 

8.28 

9.13 

9.96 

10.47 

10.55 

10.48 

10.47 

10.44 

10.09 

9.33 

8.46 

8.00 

a31 

8.21 

9.07 

9.93 

10.47 

10.57 

10.52 

10.50 

10.45 

10.07 

9.28 

8.39 

7.93 

a35 

8.29 

9.13 

9.96 

10.47 

10.55 

10.47 

10.46 

10.43 

10.09 

9.34 

8.46 

8.01 

blO 

8.27 

9.12 

9.95 

10.47 

10.55 

10.48 

10.47 

10.44 

10.09 

9.33 

8.45 

7.99 

93 

8.25 

9.10 

9.94 

10.47 

10.56 

10.50 

10.48 

10.44 

10.08 

9.31 

8.42 

7.97 

h8 

8.27 

9.12 

9.95 

10.47 

10.55 

10.48 

10.47 

10.44 

10.09 

9.32 

8.45 

7.99 

14 

8.29 

9.14 

9.96 

10.47 

10.54 

10.47 

10.46 

10.43 

10.09 

9.34 

8.46 

8.01 

18 

8.31 

9.15 

9.97 

10.47 

10.54 

10.46 

10.45 

10.43 

10.10 

9.35 

8.48 

8.03 

m18 

8.30 

9.15 

9.96 

10.47 

10.54 

10.47 

10.46 

10.43 

10.10 

9.35 

8.48 

8.03 

m23 

8.30 

9.14 

9.96 

10.47 

10.54 

10.47 

10.46 

10.43 

10.10 

9.34 

8.47 

8.02 

m6 

8.34 

9.17 

9.98 

10.46 

10.53 

10.45 

10.44 

10.43 

10.11 

9.37 

8.51 

8.07 

n15 

8.36 

9.18 

9.98 

10.46 

10.52 

10.44 

10.43 

10.42 

10.11 

9.38 

8.53 

8.08 

s27 

8.31 

9.15 

9.97 

10.47 

10.54 

10.46 

10.45 

10.43 

10.10 

9.35 

8.48 

8.03 

s5 

8.31 

9.15 

9.96 

10.47 

10.54 

10.47 

10.46 

10.43 

10.10 

9.35 

8.48 

8.03 

t6 

8.32 

9.16 

9.97 

10.47 

10.53 

10.46 

10.45 

10.43 

10.10 

9.36 

8.49 

8.05 

ul 1 

8.36 

9.18 

9.98 

10.46 

10.52 

10.44 

10.43 

10.42 

10.11 

9.38 

8.53 

8.08 

u6 

8.33 

9.17 

9.98 

10.47 

10.53 

10.45 

10.44 

10.43 

10.10 

9.36 

8.51 

8.06 

v13 

8.31 

9.15 

9.96 

10.47 

10.54 

10.46 

10.45 

10.43 

10.10 

9.35 

8.48 

8.03 

v6 

8.34 

9.17 

9.98 

10.46 

10.53 

10.45 

10.44 

10.43 

10.11 

9.37 

8.51 

8.07 

z2 

8.31 

9.15 

9.97 

10.47 

10.54 

10.46 

10.45 

10.43 

10.10 

9.35 

8.48 

8.03 

z4 

8.26 

9.11 

9.95 

10.47 

10.56 

10.49 

10.48 

10.44 

10.09 

9.32 

8.43 

7.98 
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Paso 3. 
Modelo de Angstrom 


El modelo parte del supuesto ffsico que hay una relacion directa entre la radiacion solar global y el 
brillo solar, llegando a establecer: 

Y = a + bx Ec. 3-14 

En donde: 

H 

Y = — (Radiacion medida entre Radiacion extraterrestre) 

H 0 

Hi 

x = — (Duracion medida del dfa entre duracion astronomica del dfa) 

N d 

a y b: son los coeficientes de regresion, llamados coeficientes de Angstrom. 


Para una estacion determinada se tienen los valores H y Ni (promedios mensuales obtenidos 
mediante mediciones historicas de campo). N d se calcula mediante la ecuacion 3-10 y H 0 se calcula 
mediante la ecuacion 3-13. Con todos los pares ordenados H, H 0 , Ni y N d para los diferentes meses, 
para cada estacion meteorologica se realiza un ajuste lineal para encontrar los coeficientes ay b. La 
ecuacion a resolver por mfnimos cuadrados es: 



Ec. 3-15 


Por ejemplo, para la estacion de Ahuachapan calculamos los pares ordenados (X, Y) para cada uno 
de los meses del ano, ver tabla 3-4, posteriormente hacemos una regresion lineal de los datos para 
encontrar las constantes a y b, como se muestra en la figura 3-2, los valores determinados son: b= 
0.418 y a = 0.252. 


Tabla 3-4 Calculo de constantes de Angstrom para la estacion H-8. 3 



Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

H 

5.1 

5.5 

5.8 

5.6 

5.2 

4.9 

5.6 

5.3 

4.6 

4.8 

4.9 

4.8 

N 

9.5 

9.7 

9.6 

9.0 

7.5 

6.7 

8.3 

7.4 

5.9 

7.1 

8.6 

9.4 

Ho 

8.27 

9.12 

9.95 

10.47 

10.55 

10.48 

10.47 

10.44 

10.09 

9.32 

8.45 

7.99 

Nd 

11.27 

11.56 

11.92 

12.31 

12.65 

12.81 

12.74 

12.45 

12.07 

11.68 

11.35 

11.19 

Y 

0.6139 

0.6009 

0.5843 

0.5365 

0.4904 

0.4718 

0.5309 

0.5079 

0.4570 

0.5138 

0.5832 

0.6046 

X 

0.8426 

0.8388 

0.8054 

0.7308 

0.5931 

0.5230 

0.6517 

0.5942 

0.4887 

0.6079 

0.7578 

0.8399 


Q 

El procedimiento de calculo ilustrado en la tabla 3-4 es correcto; sin embargo, las estimaciones para H8 no 
incluyen el estado actual de la base de datos digitalizada de radiacion solar y brillo solar. Los valores correctos 
de a y b, con la informacion actualizada, para la estacion H8 se presentan en la tabla 3-5. 
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Figura 3-2 Regresion por mfnimos cuadrados para los datos de Ahuachapan. 


El procedimiento se repite para las nueve estaciones que tenemos mediciones de radiacion solar y 
brillo solar, obteniendo los valores mostrados en la tabla 3-5. Para comprobar los calculos se 
sumaron las constantes “a” y “b” y estos valores oscilan alrededor de 0.65 como recomienda la 
literatura para estas latitudes y climas tropicales como el nuestro. 


Tabla 3-5 Constantes “a” y “b” de Angstrom y 
su coeficiente de correlacion “r” 


Estacion 

b 

a 

’a+b’ 

r 

a31 

0.524 

0.182 

0.706 

0.942 

93 

0.351 

0.303 

0.654 

0.952 

h8 

0.445 

0.213 

0.658 

0.966 

n15 

0.304 

0.290 

0.594 

0.981 

s27 

0.449 

0.266 

0.715 

0.973 

s5 

0.351 

0.308 

0.659 

0.964 

u 1 1 

0.359 

0.273 

0.632 

0.959 

v13 

0.309 

0.287 

0.596 

0.985 

z4 

0.567 

0.147 

0.714 

0.938 


Paso 4. 


Evaluacion de la energi'a solar disponible para lugares 
que poseen o solo brillo solar o solo radiacion solar global. 

De las 22 estaciones meteorologicas 12 poseen registro de brillo solar unicamente y 1 tiene solo 
registro de radiacion solar global. Para evaluar el potencial solar en estos casos se asignaron a estos 
lugares las constantes angstrom de las otras estaciones que si poseian registros de radiacion y brillo. 
El criterio utilizado se baso en los siguientes aspectos: 

■ Zona de vida en que estaban situadas las estaciones con informacion incompleta. 

■ La elevacion sobre el nivel del mar. 
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La cartografia de zonas de vida de El Salvador 4 describe zonas ecologicas basadas en las 
elevaciones, microclimas y temperaturas promedios, registrandose en nuestro pais 10 zonas de vida, 
as! como se muestra en la tabla 3-6 y figura 3-3. 




Mind mi dt 
V ■ 




swera 


Si las estaciones sin registros de radiacion se encuentran en la misma zona de vida, ver figura 3-4, 
que otra estacion cuya disponibilidad de energia solar se conoce y sus elevaciones sobre el nivel del 
mar son comparables, se aplicaran los coeficientes a y £> conocidos. Los resultados se muestran en la 
tabla 3-7. 



□ bh-s □ bmh-MBS [J bs-T J bmh-S □ bmh-MS □ bh-S(c) 

□ bh-MBS BAgua □ bh-T | bh-MB Qbh-T(c) 


Figura 3-3 Mapa de Zonas de Vida de El Salvador. Puntos rojos muestran estaciones con registros 

de radiacion y los verdes con registro de brillo. 


Tabla 3-6 Zonas de Vida de El Salvador 


Simbolo 

Descripcion 

bh-MB 

Bosque humedo Montano Bajo 

bh-MBS 

Bosque humedo Montano Bajo Subtropical 

bh-S 

Bosque humedo Subtropical 

bh-S(c) 

Bosque humedo Subtropical 

bh-T 

Bosque humedo Tropical, transicion a Subtropical 

bmh-MBS 

Bosque muy humedo Montano Bajo Subtropical, transicion a humedo 

bmh-MS 

Bosque muy humedo Montano Subtropical 

bmh-S 

Bosque muy humedo Subtropical, transicion a humedo 

bs-T 

Bosque seco Tropical 

bs-T(c) 

Bosque seco Tropical, transicion a Subtropical (con biotemperatura < 
24 2 C.) 


4 La zona de vida elaborado por L. R. Holdridge y ampliado por Joseph Tosi y Gary Hartshorn, CATIE y la 
Direccion General de Recursos Naturales Renovables del Ministerio de Agricultura y Ganaderfa. Este mapa se 
encuentra disponible en el Sistema Ambiental de El Salvador, publicado por el MARN. 
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Figura 3-4 Ubicacion de las estaciones con respecto a las zonas de vida. Los puntos en amarillo son 
aquellas estaciones con coeficientes a y b segun el modelo de Angstrom, los verdes los que no lo 

poseen y las flechas denotan la asignacion realizada. 


Tabla 3-7 Coeficientes del modelo de Angstrom. Datos del tipo 1 ya poseen registros de radiacion y 
los del tipo 2 se les asignado las a y b de las estaciones mostradas en la columna fuente. 


Estacion 

Tipo 

Elevacion 

msnm 

Fuente 

H-8 

1 

725 


A-18 

2 

1770 

A-31 

A-31 

1 

1851 


A-35 

2 

2027 

Z-4 

T-6 

2 

15 

U-11 

L-4 

2 

460 

V-13 

L-8 

2 

960 

H-8 

S-5 

1 

700 


S-27 

1 

650 


G-3 

1 

320 


G-13 

1 

1960 

(A-31) 

V-6 

2 

30 

U-11 

V-13 

1 

570 


U-6 

2 

920 

Z-4 

U-11 

1 

75 


M-6 

2 

80 

N-15 

M-18 

2 

195 

N-15 

M-23 

2 

1370 

Z-4 

Z-2 

2 

250 

N-15 

Z-4 

1 

1900 


N-15 

1 

95 


B-10 

2 

245 

N-15 
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N d y H 0 se calculan en la forma ya explicada anteriormente y a y b se obtienen de acuerdo a los 

criterios anteriores, de manera que la radiacion solar disponible sera la calculada por la ecuacion 3- 
16, en el caso en que la estacion no posea registro de brillo (como la G-13) la ecuacion a utilizar es la 
3-17. 


H = a + b 


f N: ^ 


\ N d J 


H 


o 


Ec. 3-16 



f H-a' 


V 


H 


o 


J 



Ec. 3-17 


En donde los diferentes terminos fueron descritos anteriormente. 


3.3 Datos de radiacion v brillo solar medidos v calculados por el modelo 

de Angstrom para estimation de los mapas de potencial solar. 


Tabla 3-8 Datos utilizados para estimacion del mapa de radiacion solar (KWh/m 2 /dia) (Datos de 

Solarimetria + SWERA) . 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a18 

4.2 

4.9 

5.4 

5.3 

4.3 

3.9 

4.8 

4.6 

3.7 

3.6 

3.7 

3.6 

a31 

4.9 

5.2 

5.0 

4.5 

3.6 

4.1 

5.2 

4.7 

3.8 

4.4 

4.6 

4.6 

a35 

4.4 

4.7 

4.7 

4.4 

4.0 

3.7 

5.1 

5.0 

3.7 

3.7 

4.2 

4.2 

blO 

4.4 

4.8 

5.3 

5.2 

4.9 

4.7 

5.1 

5.1 

4.5 

4.5 

4.3 

4.2 

g13 

4.6 

5.0 

5.6 

5.3 

5.1 

5.0 

5.5 

5.5 

4.8 

4.7 

4.3 

4.3 

q3 

4.9 

5.4 

5.7 

5.5 

5.2 

5.2 

5.4 

5.4 

5.0 

4.9 

4.7 

4.6 

h8 

5.0 

5.4 

5.7 

5.5 

5.0 

4.7 

5.2 

5.0 

4.3 

4.5 

4.6 

4.5 

14 

4.5 

4.9 

5.1 

5.1 

4.8 

4.6 

5.1 

4.9 

4.4 

4.3 

4.3 

4.3 

18 

5.0 

5.4 

5.7 

5.3 

4.7 

4.4 

5.2 

5.0 

4.4 

4.4 

4.7 

4.7 

m18 

4.6 

5.1 

5.4 

5.3 

4.9 

4.8 

5.1 

5.0 

4.6 

4.4 

4.4 

4.3 

m23 

4.6 

5.2 

5.4 

4.5 

3.7 

3.7 

4.9 

4.7 

4.0 

3.9 

4.2 

4.0 

m6 

4.5 

5.0 

5.4 

5.4 

5.0 

4.9 

5.2 

5.2 

4.8 

4.6 

4.4 

4.4 

n15 

4.7 

5.1 

5.4 

5.3 

4.8 

5.0 

5.2 

5.2 

4.8 

4.7 

4.5 

4.4 

s27 

5.1 

5.6 

6.2 

5.9 

5.2 

5.3 

5.9 

5.6 

4.9 

4.8 

5.0 

4.8 

s5 

4.9 

5.4 

5.7 

5.4 

4.9 

5.1 

5.5 

5.2 

4.6 

4.8 

4.8 

4.8 

t6 

4.9 

5.3 

5.6 

5.5 

5.1 

4.9 

5.3 

5.2 

4.8 

4.7 

4.7 

4.6 

ul 1 

4.9 

5.4 

5.5 

5.2 

4.8 

4.8 

5.4 

5.4 

5.0 

4.8 

4.8 

4.7 

u6 

5.0 

5.6 

5.8 

5.4 

4.6 

4.2 

5.2 

5.1 

4.1 

4.2 

4.6 

4.8 

v13 

4.6 

5.2 

5.3 

5.1 

4.8 

4.7 

5.1 

5.1 

4.6 

4.5 

4.6 

4.4 

v6 

4.8 

5.3 

5.6 

5.5 

5.1 

4.9 

5.3 

5.3 

4.8 

4.7 

4.6 

4.6 

z2 

4.5 

5.0 

5.3 

5.2 

4.8 

4.8 

5.1 

5.1 

4.7 

4.6 

4.4 

4.3 

z4 

4.4 

4.9 

5.1 

4.8 

4.1 

4.0 

4.5 

4.3 

3.5 

3.7 

3.7 

3.9 
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Tabla 3-9 Datos utilizados para estimacion del mapa de brillo solar (horas de luz solar) (Datos 

Solarimetna + SWERA). 


Estacion 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

a18 

6.9 

7.8 

8.2 

7.6 

5.5 

4.7 

6.7 

6.0 

4.3 

4.6 

5.4 

5.8 

a31 

8.7 

8.5 

8.1 

7.1 

5.0 

4.8 

6.7 

6.2 

4.0 

5.9 

7.7 

8.0 

a35 

7.7 

7.5 

6.9 

5.9 

5.1 

4.7 

7.6 

7.2 

4.6 

5.2 

7.0 

7.5 

blO 

8.9 

9.0 

9.4 

8.5 

7.0 

6.7 

8.1 

8.0 

6.2 

7.3 

8.3 

8.7 

g13 

8.2 

8.0 

8.6 

7.6 

7.1 

7.2 

8.4 

8.2 

6.8 

7.1 

7.1 

7.6 

93 

9.1 

9.6 

9.2 

8.5 

7.5 

6.8 

7.9 

7.4 

6.3 

7.3 

8.2 

8.1 

h8 

9.1 

9.4 

9.4 

9.0 

7.4 

6.6 

8.2 

7.4 

5.9 

7.1 

8.6 

9.4 

14 

9.2 

9.2 

8.8 

8.0 

6.8 

6.4 

8.1 

7.4 

5.7 

6.7 

8.3 

8.9 

18 

9.8 

9.8 

9.5 

8.1 

6.7 

6.0 

8.0 

7.4 

5.9 

6.6 

8.6 

9.5 

m18 

9.6 

10.0 

9.7 

8.6 

7.3 

7.0 

8.3 

7.6 

6.3 

7.0 

8.4 

9.2 

m23 

8.0 

8.6 

8.4 

6.2 

4.5 

4.7 

7.2 

6.6 

5.3 

5.5 

6.9 

7.0 

m6 

9.4 

9.8 

9.7 

8.9 

7.7 

7.3 

8.6 

8.6 

7.2 

7.8 

8.5 

9.3 

n15 

9.8 

10.0 

9.7 

8.9 

7.4 

7.4 

8.8 

8.5 

7.2 

8.1 

9.1 

9.7 

s27 

8.5 

9.1 

9.2 

8.5 

6.4 

6.4 

8.2 

7.7 

5.9 

6.9 

7.9 

8.2 

s5 

9.1 

8.7 

9.2 

7.5 

5.6 

5.5 

7.6 

7.3 

5.6 

7.3 

8.3 

9.0 

t6 

9.8 

9.9 

9.6 

8.6 

7.3 

6.8 

8.2 

7.7 

6.7 

7.3 

8.9 

9.4 

ul 1 

9.7 

9.9 

9.7 

8.6 

7.0 

6.9 

8.2 

8.2 

7.0 

7.3 

8.9 

9.5 

u6 

9.1 

9.4 

9.1 

7.9 

6.5 

5.6 

7.8 

7.5 

5.5 

6.2 

7.9 

8.9 

v13 

9.8 

10.0 

9.5 

8.5 

6.9 

6.5 

8.1 

8.0 

6.6 

7.2 

8.9 

9.7 

v6 

9.4 

9.7 

9.5 

8.5 

7.4 

6.9 

8.3 

8.2 

6.8 

7.4 

8.5 

9.4 

z2 

9.4 

9.8 

9.4 

8.5 

7.0 

7.0 

8.2 

8.0 

6.9 

7.6 

8.6 

9.2 

z4 

7.0 

7.9 

7.9 

7.1 

5.6 

5.1 

7.1 

6.2 

4.6 

4.7 

5.7 

6.2 
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4. Metodo de Interpolacion de Kriging para estimacion 

del potencial solar en El Salvador. 

Una vez aplicado el modelo de Angstrom a las estaciones meteorologicas del pais se obtuvo un total de 
22 estaciones con datos calculados y estimados de radiacion solar y brillo solar para cada uno de los 
meses del ano y un promedio anual para cada uno de estos dos parametros. Los valores de radiacion y 
brillo solar se muestran en las tablas 3-8 y 3-9 respectivamente. 

Luego se evaluaron diferentes metodos de interpolacion con los datos obtenidos descartando el Inverse 
Distance Weighted (IDW) principalmente por el patron de forma “duck-egg” (huevo de pato) que genera 
alrededor de puntos de datos solitarios cuyos valores difieren grandemente del de sus vecinos. Tambien 
este metodo de interpolacion (IDW) trabaja mejor cuando los datos muestreados son numerosos y estan 
distribuidos uniformemente lo cual no era nuestro caso. Ademas se analizo el metodo de interpolacion de 
Spline y se descarto por la siguiente razon: el metodo asume que una funcion aproximada debe pasar tan 
cerca como sea posible de los puntos de datos y a su vez sea lo mas suave posible; esto produce que 
algunos puntos estimados queden muy por arriba o por debajo de la funcion suave estimada. En ambos 
casos los pesos asignados para calcular la combinacion lineal de la aproximacion no son optimos. Debido 
a estas razones y dado que en el metodo de interpolacion de Kriging asigna pesos de manera optima se 
opto utilizarlo para el caso del mapa de brillo solar y radiacion solar. 


4.1 Metodo de interpolacion de Kriging. 


Para propositos academicos se describe a continuacion el metodo de interpolacion de Kriging. Cuando se 
comparan los metodos de interpolacion IDW y Spline se observa que retornan modelos con grandes 
diferencias en los patrones globales, en la cantidad de detalle local y en las predicciones de los valores 
maximos y rm'nimos interpolados. En estos metodos, la unica manera de determinar “lo bueno” de las 
predicciones es calcular estimados para un conjunto extra de puntos de validation que no hayan sido 
usados en la interpolacion original. Otra desventaja de estos metodos es que no hay un metodo a priori 
para saber si se han seleccionado los mejores valores para los pesos o si el tamano para la busqueda de 
los vecinos es la apropiada. Estas desventajas motivaron al geomatematico Frances Georges Matheron y 
al ingeniero de minas sudafricano D.G. Krige a desarrollar metodos optimos de interpolacion para uso en 
la industria de minas. Estos metodos de interpolacion optima basados en Geoestadistica se usan 
ampliamente para modelar agua subterranea y mapeo de suelos. 

Los metodos Geoestadisticos para interpolacion comienzan con el reconocimiento de que la variacion 
espacial de cualquier atributo continuo es a menudo demasiado irregular para ser modelada por una 
simple funcion matematica suave. En lugar de esto, la variacion puede describirse de mejor manera por 
medio de una superficie estocastica. El atributo es conocido entonces como una variable 
regionalizada; el termino se aplica por igual a la variacion de la presion atmosferica, a la elevation sobre 
el nivel del mar o a la distribution de indicadores continuos demograficos. La interpolacion con 
Geoestadistica es conocida como Kriging en honor de D.G. Krige. Este metodo es optimo en el 
sentido de que los pesos de la interpolacion se seleccionan de tal forma que optimicen la funcion de 
interpolacion. 
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Teona de variable regionalizada . 1 

La teoria de variable regionalizada asume que la variacion espacial de cualquier variable puede 
expresarse como la suma de tres componentes principales. Estos son: (a) un componente estructural, 
que tiene una media o tendencia constante; (b) Un componente aleatorio, pero espacialmente 
correlacionado conocido como la variacion de la variable regionalizada, y (c) ruido aleatorio no 
correlacionado espacialmente tambien llamado termino de error residual. 





( 3 ) 


x 


Figura 4-1 Teoria de variable regionalizada: (1) diferencias en niveles promedio (a) y tendencia o drift (b); 
(2) espacialmente correlacionada pero variacion irregular (aleatoria); (3) variacion local no correlacionada 

y aleatoria. 

Sea x una posicion en tres dimensiones. Entonces el valor de una variable aleatoria Z en x se expresa 
como: 


Z (x) = m(x) + s' (x) + s" (A) 

en donde m(x) es la funcion determimstica que describe el componente “estructural” de Z en x, s’( x ) es el 
termino que denota el componente estocastico de m(x) - la variable regionalizada y s” es un ruido 
Gaussiano independiente espacialmente que tiene media cero y varianza o 2 . Note que el uso de la 
mayuscula Z indica que es una funcion aleatoria y no un atributo medido z. 

El primer paso es decidir en una funcion apropiada para m(x). En el caso mas simple, en donde no existe 
tendencia o “drift”, m(x) es igual al valor medio en el area de muestreo y el promedio o diferencia 
esperada entre cualquiera dos lugares x y x + h separados por una distancia vectorial h, sera cero: 

E[Z(x)-Z(x + h)\ = 0 

en donde Z(x), Z(x+h) son los valores de variable aleatoria Z en los lugares x y x+h. Tambien se asume 
que la varianza de las diferencias depende unicamente de la distancia entre los sitios, h, tal que: 


1 La teoria de variable regionalizada y las graficas de esta seccion han sido adaptadas del libro “Principles of 

Geographical Information Systems” de Peter Burrough y Rachael McDonell, Oxford University Press, 1998. 
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{Z(jc) - Z(jc + h)f = E {s'ix) -e\x + h)} 2 = 2 y(h) 


En donde y (h) es conocida como la semivarianza. 

Las dos condiciones, estacionaridad de diferencias y varianza de diferencias, definen los requerimientos 
para la hipotesis intrinseca de la teoria de variable regionalizada. Esto significa que una vez los efectos 
estructurales han sido tornados en cuenta, la variacion restante es homogenea en su variacion de tal 
forma que las diferencias entre sitios son sencillamente una funcion de la distancia entre ellos. 
Reescribimos la ecuacion (A) como: 


Z(x) = m(x) + y (h) + s" 

Para mostrar la equivalencia entre s' (x) y y ( h ) se hace la consideracion siguiente: Si las condiciones 

especificadas por la hipotesis intrinseca se cumplen, la semivarianza puede ser estimada de los datos 
muestreados por: 


/(h) = 2 ^^ {z (Xj) - z ( Xj+h )} 2 Ecuacion 4.1 .A 

i = 1 

en donde n es el numero de pares de puntos de observaciones de los valores del atributo z separados 
por una distancia h. 

Al grafico de y(h) en funcion de h es conocido como el “variograma experimental”. 

El variograma proporciona informacion util para la interpolacion, muestreo optimizado y para determinar 
patrones espaciales. Para hacer esto, sin embargo, debemos primero ajustar un modelo teorico al 
variograma experimental. 


Ajustando modelos de variogramas 

La figura 4-2 muestra un variograma experimental tipico. Para valores grandes de “lag”, h, la curva se 
nivela. Esta parte horizontal es conocida como el “sill”; implica que en estos valores del lag no hay 
dependencia espacial entre los puntos de datos, debido a que todos los estimados de las varianzas de 
las diferencias seran invariantes con la distancia de separacion de las muestras. 
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La curva crece de un valor bajo de y (h) hasta el sill, alcanzandolo a un valor de h conocido como el 
rango. Esta es la parte critica e importante del variograma debido a que describe como las diferencias 
entre sitios son dependientes espacialmente. El rango responde a la pregunta de que tan grande 
debe ser la ventana de interpolacion. Claramente, si la distancia que separa a un sitio no visitado 
de un punto de datos es mayor que el rango, entonces el punto de datos no hace una contribucion 
util a la interpolacion. El tercer punto mostrado en la figura es el que en el modelo ajustado no pasa a 
traves del origen, sino que corta el eje Y en un valor positivo de y (h) cuando h tiende a 0 es un 
estimado de e", el residual o ruido espacialmente no correlacionado. e" es conocido como el “nugget”. 

La forma del variograma es bastante reveladora acerca de la clase de variacion espacial presente en un 
area, y puede ayudar a decidir como proceder con la interpolacion. Cuando la varianza del nugget es 
considerable pero no demasiado grande y hay claridad en el rango y en el sill, una curva conocida como 
el modelo esferico se podria ajustar al variograma observado. 


Modelo esferico 


Y (h) = Co 


r 



2a 


V 


1 fh' 

2\a y 



> para 0 < h < a 


= C 0 + Ci para h > a 
Y (0) = 0 


en donde “a” es el rango, h es el lag, C 0 es la varianza del nugget y C 0 + Ci es igual al sill. 

Si hay un claro nugget y sill, pero unicamente una aproximacion gradual al rango, el modelo exponencial 
es la mejor seleccion: 


y (h) - Co + Ci 


1-exp 




Si la variacion es muy suave y la varianza del nugget e" es muy pequena comparada con la variable 
aleatoria e" (x), entonces el variograma se ajusta mejor a una curva que tiene una inflexion tal como el 
modelo Gaussiano. 


y (h) - Co + Ci 



Todos estos modelos son conocidos como varioqramas transitivos debido a que la correlacion espacial 
varia con h. Los variogramas no transitivos no tienen sill dentro del area muestreada y podrian ajustarse 
al modelo lineal. 


y ( h ) = C 0 + bh 


Ecuacion 4.1. B 


en donde b es la pendiente de la recta. 
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Figura 4-3 Ejemplos de modelos de variogramas. 


Usando el variograma para interpolation : Kriging Ordinario 

El variograma ajustado se puede usar para determinar los pesos ^ i necesitados para la interpolation 
local en un punto desconocido Xo. 

n n 

Z (x 0 ) = V A i z(Xi) con V/li = 1 Ecuacion 4.1. C 

i = 1 i = 1 


Los pesos'^ se seleccionan de tal forma que el estimado Z (x 0 ) sea no sesgado (“ unbiased”) y que la 

— 2 

varianza de la estimation <j e sea menor que cualquier combination lineal de los valores observados. 
La minima varianza de Z (x D ), esta dada por: 


G e = y{ Xi , x 0 ) +</> 

i = 1 

y se obtiene cuando : 


Ecuacion 4.1. D 


2>r(Xi.XiW = ^(Xj.Xo) Vj Ecuacion 4.1. E 

1=1 

La cantidad y ( Xj , Xj ) es la semi-varianza de z entre los puntos de datos Xj y Xj ; y ( Xi , x 0 ) es la semi- 
varianza entre el punto de datos Xj y el punto no visitado x 0 Ambas cantidades se obtienen del 
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variograma ajustado. La cantidad (/) es un multiplicador de Lagrange requerido para la minimizacion. 

Este metodo es conocido como Kriqinq Ordinario . Es un interpolador exacto en el sentido que los 
valores interpolados coinciden con los valores en los puntos de datos. 


— 2 

La estimacion del error <j e conocido como varianza de Kriging, tambien se puede mapear y proporcionar 
informacion a cerca de la confiabilidad de los valores interpolados en el area de interes. A menudo la 
varianza de Kriging se mapea como desviacion estandar Kriging (o error Kriging) debido a que tiene las 
mismas unidades que las predicciones. 

Ejemplo de interpolacion de Kriging. 


Procedimiento de Kriging en terminos matriciales. 

La ecuacion 4.1.E pueden escribirse en terminos matriciales como sigue: 

V*A = g 


Yu Yin V 





• • • • 

• • • • 

Ym Y nn 1 

1 ... 1 0 y 

* 

• 

• 

1<P J 

— 

• 

• 

Yn 0 

l 1 J 


La matriz gamma T contiene los valores del semivariograma modelado entre todos los pares de 
localizaciones con valores conocidos. Los y tj denotan los valores del semivariograma modelado 

calculados en base a la distancia entre las localizaciones i y j. El vector g contiene los valores del 
semivariograma modelado entre cada localizacion conocida la localizacion a predecir (N indica el numero 
de localizaciones con valores conocidos y 0 es la localizacion a estimar su valor). Durante los calculos 
matriciales se estima el valor de cp que surge como resultado de la restriccion de que los pesos son no 

n 

sesgados ( ^ i A\ = 1 ). El objetivo de la interpolacion de Kriging es resolver la ecuacion A=r~ 1 *g 

i = 1 

para luego estimar el valor en la localizacion a predecir, usando la ecuacion 4.1 .C. 


Ejemplo numerico. 

En la tabla 4-1 aparecen las estaciones con valores conocidos (radiacion promedio anual del mes de 
enero) y su ubicacion geografica en el pais. Se desea estimar el valor en el punto descrito en la tabla 4-2. 


Tabla 4-1 Cuadro de Coordenadas y Valores conocidos. 


Estacion 

Coordenada 

X 

Coordenada 

Y 

Valor de 
Radiacion 
de Enero 

m18 

568868.255 

288152.772 

4.6 

m23 

583095.785 

293725.708 

4.6 

m6 

594472.182 

257998.802 

4.5 

u6 

557209.568 

262868.473 

5.0 

z2 

596640.804 

286393.651 

4.5 
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Tabla 4-2 Coordenada del punto de Interes a estimar por el metodo de Kriging 


Punto de Interes 

X 

Y 

A 

585548 

280009 


Paso 1 . El primer paso para crear el semivariograma empirico es calcular la distancia y la diferencia al 
cuadrado entre cada par de localizaciones conocidas. La distancia entre dos localizaciones se calcula 

por la formula a. = J(x ; -x ) + (y,— y.) . La semivarianza empirica se calcula por la ecuacion 

4.1. A, es decir, 0.5 * promedio(valor_en_i — valor _en_j) . Los resultados se presentan en la 
tabla 4-3. 


Tabla 4-3 Cuadro de distancias y calculos del paso 1 . 


Estl 

Est2 

Distancia (m) 

Diferencias 
al cuadrado 

Semi- 

varianza 

m18 

m23 

15280.06 

0 

0 

m18 

m6 

39557.84 

0.01 

0.005 

m18 

u6 

27842.79 

0.16 

0.08 

m18 

z2 

27828.20 

0.01 

0.005 

m23 

m6 

37494.46 

0.01 

0.005 

m23 

u6 

40277.35 

0.16 

0.08 

m23 

z2 

15402.16 

0.01 

0.005 

m6 

u6 

37579.46 

0.25 

0.125 

m6 

z2 

28477.54 

0 

0 

u6 

z2 

45915.75 

0.25 

0.125 


Paso 2. Usando un tamafio de lag de 9000 metros se ubican las distancias y los pares de puntos 
conocidos dentro de cada lag. Se estima el promedio de las distancias entre pares y el promedio de las 
semivarianzas de la tabla 4-3. Los resultados se presentan en la tabla 4-4. 


Tabla 4-4 Cuadro de Lags y calculo se semivariograma empirico (tecnica de Binning). 


# 

Lags 

Distancia (m) 

Vectores en el rango de distancia 

Distancia 
promedio de 
los pares 

Semi- 

varianza 

promedio 

Desde 

Hasta 

1 

9000 

18000 

(ml 8,m23),(m23,z2) 

15341.11 

0.0025 

2 

18000 

27000 




3 

27000 

36000 

(m18,u6),(m18,z2),(m6,z2) 

28049.51 

0.02833 

4 

36000 

45000 

(m18,m6),(m23,m6),(m23,u6),(m6,u6) 

38727.28 

0.05375 

5 

45000 

54000 

(u6,z2) 

45915.75 

0.125 


Paso 3. Se grafica los pares (distancia promedio de pares, semivarianza promedio) y se ajusta un modelo 
lineal de semivariograma (ver ecuacion 4.1.B). El modelo a utilizar es 0.000002*/?, en donde h es la 
distancia entre pares de puntos. 
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Figura 4-3-a Grafico de Semivariograma empirico y modelo ajustado. 


Paso 4. Calculo de la matriz gamma. Para ello primero se calcula la matriz de distancias y se evalua cada 
elemento de la matriz con el modelo de semivariograma ajustado, en nuestro caso 0.000002* h . 

Tabla 4-5 Matriz de Distancias. 



m18 

m23 

m6 

u6 

z2 

m18 

0.0 

15280.1 

39557.8 

27842.8 

27828.2 

m23 

15280.1 

0.0 

37494.5 

40277.4 

15402.2 

m6 

39557.8 

37494.5 

0.0 

37579.5 

28477.5 

u6 

27842.8 

40277.4 

37579.5 

0.0 

45915.8 

z2 

27828.2 

15402.2 

28477.5 

45915.8 

0.0 


Tabla 4-6 Construccion de Matriz Gamma. 



m18 

m23 

m6 

u6 

z2 

Lagrange 

m18 

0.0000 

0.0306 

0.0791 

0.0557 

0.0557 

1 

m23 

0.0306 

0.0000 

0.0750 

0.0806 

0.0308 

1 

m6 

0.0791 

0.0750 

0.0000 

0.0752 

0.0570 

1 

u6 

0.0557 

0.0806 

0.0752 

0.0000 

0.0918 

1 

z2 

0.0557 

0.0308 

0.0570 

0.0918 

0.0000 

1 

Lagrange 

1 .0000 

1.0000 

1.0000 

1.0000 

1.0000 

0 


Paso 5. Se encuentran los pesos X usando la ecuacion A = T 1 * g . 


Tabla 4-7 Matriz Inversa de Gamma. 



m18 

m23 

m6 

u6 

z2 

Lagrange 

m18 

-23.7931 

15.9513 

-0.1574 

8.3011 

-0.3019 

0.0795316 

m23 

15.9513 

-30.0111 

-0.8328 

-0.7631 

15.6556 

0.1541849 

m6 

-0.1574 

-0.8328 

-11.8958 

5.0383 

7.8476 

0.2492623 

u6 

8.3011 

-0.7631 

5.0383 

-11.4389 

-1.1374 

0.3249942 

z2 

-0.3019 

15.6556 

7.8476 

-1.1374 

-22.0640 

0.1920267 

Lagrange 

0.0795 

0.1542 

0.2493 

0.3250 

0.1920 

-0.0532174 
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Tabla 4-8 Calculo del vector g. 


Punto 

Distancia con 
el punto de 
interes 

gVector 

m18 

18561.65 

0.0371233 

m23 

13934.18 

0.02786836 

m6 

23750.58 

0.04750115 

u6 

33118.94 

0.06623789 

z2 

12798.99 

0.02559797 



1 


Tabla 4-9 Vector de pesos A : Inversa de Gamma * Vector g. 


Estacion 

Pesos ( A ) 

m18 

0.1754 

m23 

0.2206 

m6 

0.1898 

u6 

0.0644 

z2 

0.3498 

(p 

-0.008 


Paso 6. Se calcula el valor a estimar usando la ecuacion 4.1.C y el error estandar de Kriging (tablas 4-10 
y 4-11). 


Tabla 4-10 Prediccion de Kriging para el punto A. 


Estacion 

Radiacion 
del mes de 
Enero 

Aporte de 
cada 
estacion 

m18 

4.6 

0.807 

m23 

4.6 

1.015 

m6 

4.5 

0.854 

u6 

5.0 

0.322 

z2 

4.5 

1.574 

Valor estimado 

4.572 


Tabla 4-1 1 Varianza y error estandar de Kriging. 


Pesos ( A ) 

Vector g 

g*A 

0.17543 

0.03712 

0.00651 

0.22064 

0.02787 

0.00615 

0.18976 

0.04750 

0.00901 

0.06441 

0.06624 

0.00427 

0.34976 

0.02560 

0.00895 

-0.00769 

1.00000 

-0.00769 

Varianza de Kriging 

0.02721 

Error estandar de Kriging 

0.16495 
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Figura 4-3-b Resultados de la estimacion de Kriging. 

El punto en anaranjado con valor 4.572 es el valor estimado en esa localizacion. En cada uno de los 
puntos conocidos se indican los pesos con los que contribuyen a la estimacion: ml 8 (0.1754), m23 
(0.2206), m6 (0.1898), u6 (0.0644) y z2 (0.3498). De la figura 4-3-b es claro que el punto que mas 
contribuye a la estimacion es z2 con un valor de 4.5 y el punto que menos contribuye es u6 con un valor 
de 5.0. 


4.2 Calculo de los promedio anuales de radiacion solar v brillo solar 

usando Kriging ordinario. 

Para el calculo de los mapas de promedios anuales de radiacion y brillo solar y sus respectivos mapas 
mensuales se utilizaron los siguientes parametros en el modulo de analisis geoestadistico de ArcGis: 

■ Tamano del lag igual a 9000 metros. Para ello se utilizo la regia empirica de que el lag size * # de 
lags * 0.5 * distancia maxima entre las estaciones a interpolar, es decir lag size * 0.5 * 225 Km / 
12 * 9.375 Km. Se selecciono como numero de lags el valor de 12, el cual sugiere ArcGis por 
defecto. 

■ Como modelo de semivariograma se utilizo el tipo exponencial ya que experimentalmente arrojo 
los valores mas bajos del promedio de error estandar. 

■ Para el caso del promedio anual de radiacion solar el promedio del error estandar esta en el 
rango de 0.07 a 0.36. El valor promedio del valor del error fue de 0.3005. Si se asume que los 
errores estan normalmente distribuidos, los intervalos de prediccion del 95% se calculan por la 
formula valor estimado por Kriging ± 1.96 * error estandar de Kriging. Por ejemplo, si el valor y 
error estandar estimado por Kriging en un punto de interes son 4.45 y 0.21 respectivamente, 
entonces con un 95% de probabilidad el valor de radiacion en ese punto se encuentra en el rango 
[4.04, 4.86]. 

■ Para el caso del promedio anual de brillo solar el promedio del error estandar esta en el rango de 
0.68 a 0.86. 
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Mapa de Radiacion Solar 

Promedio Anual - El Salvador 
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Figura 4-5 (a) Mapa de radiacion (promedio anual) en KWh/m 2 /dia para El Salvador. Interpolacion de Kriging. 
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Figura 4-5 (b) Mapa del error promedio estandar en radiacion (promedio anual) en KWh/m 2 /dfa para El Salvador. 
Interpolation de Kriging. 
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Figura 4-6 (a) Mapa de brillo solar (en horas de luz solar) para El Salvador. Interpolacion de Kriging. 
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Figura 4-7 (a) Horas de luz solar mensuales (enero a junio). 
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Figura 4-7 (b) Horas de luz solar mensuales (julio a diciembre). 
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Figura 4-8 (a) Radiacion solar mensual (enero a junio). 
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Figura 4-8 (b) Radiation solar mensual (julio a diciembre). 
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Tabla 4-1 Tabla de interpolation de kriging para la radiacion solar en El Salvador. 


Period 

Sill 

Nugget 

Range 

Lag Size 

Equation 

Average 

Standard 

Error 

Anual 

0.10394 

0 

53772 

9000 

0. 1 0394*Exponential(53772)+0*Nugget 

0.3005 

Enero 

0.062812 

0 

31877 

9000 

0.06281 2*Exponential(31 877)+0*Nugget 

0.2555 

Febrero 

0.067177 

0 

32809 

9000 

0.0671 77*Exponential(32809)+0*Nugget 

0.2634 

Marzo 

0 

0.10446 

102170 

9000 

0*Exponential(1 021 70)+0. 1 0446*Nugget 

0.3387 

Abril 

0 

0.1361 

102170 

9000 

0*Exponential(1 021 70)+0.1 361 *Nugget 

0.3866 

Mayo 

0.083893 

0.14855 

106680 

9000 

0.083893*Exponential(106680)+0.14855*Nugget 

0.4719 

Junio 

0.22525 

0 

15234 

8000 

0.22525*Exponential(15234)+0*Nugget 

0.5089 

Julio 

0.048243 

0 

106680 

9000 

0.048243*Exponential(1 06680)+0. 03531 9*Nugget 

0.2667 

Agosto 

0.056818 

0.044953 

106680 

9000 

0.05681 8*Exponential(1 06680)+0.044953*Nugget 

0.2961 

Septiembre 

0.21217 

0 

47035 

9000 

0.21217*Exponential(47035)+0*Nugget 

0.442 

Octubre 

0.15824 

0 

45667 

9000 

0.15824*Exponential(45667)+0*Nugget 

0.3838 

Noviembre 

0.080906 

0.029254 

67634 

9000 

0.080906*Exponential(67634)+0.029254*Nugget 

0.3145 

Diciembre 

0.057443 

0.040771 

45592 

9000 

0. 057443*Exponential(45592)+0. 040771 *Nugget 

0.3218 


Tabla 4-2 Tabla de interpolation de kriging para el brillo solar en El Salvador. 


Period 

Sill 

Nugget 

Range 

Lag Size 

Equation 

Average 

Standard 

Error 

Anual 

0.28763 

106680 

0.34483 

9000 

0.28763*Exponential(106680)+0.34483*Nugget 

0.7602 

Enero 

0.29494 

106680 

0.48436 

9000 

0.29494*Exponential(106680)+0.48436*Nugget 

0.8604 

Febrero 

0.28792 

106680 

0.3537 

9000 

0.28792*Exponential(106680)+0.3537*Nugget 

0.767 

Marzo 

0.15358 

106680 

0.3763 

9000 

0.1 5358*Exponential(1 06680)+0.3763*Nugget 

0.724 

Abril 

0 

102170 

0.68797 

9000 

0*Exponential(1 021 70)+0.68797*Nugget 

0.8691 

Mayo 

0.83864 

16396 

0 

9000 

0.83864*Exponential(16396)+0*Nugget 

0.9793 

Junio 

0.80708 

16396 

0 

9000 

0.80708*Exponential(16396)+0*Nugget 

0.9607 

Julio 

0.10127 

106680 

0.20599 

9000 

0. 1 01 27*Exponential(1 06680)+0.20599*Nugget 

0.5455 

Agosto 

0.30653 

106680 

0.23187 

9000 

0.30653*Exponential(1 06680)+0.231 87*Nugget 

0.6787 

Septiembre 

0.90693 

74002 

0 

9000 

0.90693*Exponential(74002)+0*Nugget 

0.8172 

Octubre 

0.9235 

35635 

0 

9000 

0.9235*Exponential(35635)+0*Nugget 

0.9664 

Noviembre 

0.43481 

106680 

0.59974 

9000 

0.43481 *Exponential(1 06680)+0.59974*Nugget 

0.9811 

Diciembre 

0.49805 

106680 

0.81996 

9000 

0.49805*Exponential(1 06680)+0.81 996*Nugget 

1.119 
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5. Determinacion del potencial eolico en El Salvador. 


El viento se genera por el calentamiento desigual de la superficie de la tierra realizado por el Sol. En el 
caso de El Salvador predominan tres sistemas de vientos: vientos alisios de caracter regional, vientos de 
caracter local brisas mar-tierra y vientos de caracter local brisas montana-valles. Estos sistemas de viento 
fueron explicados en detalle en el capitulo 1 de este documento. Segun los resultados de este estudio las 
velocidades instantaneas maximas (rafagas) registradas en las estaciones meteorologicas del pais estan 
en el rango de 4.5 a 7.8 m/s. Segun la escala de velocidades de la tabla 5-1 estas velocidades maximas 
corresponden a vientos de tipo 3 y 4, es decir, moderados. En general las velocidades de vientos 
promedios anuales registrados por el proyecto SWERA-UCA oscilan entre 1.12 a 3.91 m/s 1 , lo que 
equivale a vientos de tipo 1 a 3 entre ligeros y moderados. 


Tabla 5-1 Escala de velocidades de viento. 


Velocidades del viento en 
m/s 

Tipo 

Descripcion 

0.0 -0.4 

0 

Calma 

0.4- 1.8 

1 

Ligero 

1.8 -3.6 

2 


3.6 -5.8 

3 


5.8 -8.5 

4 

Moderado 

8.5-11.0 

5 

Fresco 

11.0-14.0 

6 


14.0-17.0 

7 


17.0-21.0 

8 

Temporal 

21.0-25.0 

9 


25.0-29.0 

10 

Fuerte temporal 

29.0 - 34.0 

11 


>34.0 

12 

Huracan 


Nota: Velocidades de viento a 10 m de altura Escala Beaufort (anticuada) Viento. 


5.1 Conceptos basicos de viento. 


De los estudios de potencial eolico realizados en el pais se concluye que existe gran potencial para el 
desarrollo de sistemas pequenos de energi'a. La implementacion de granjas y/o parques eolicos requiere 
de mediciones mas precisas. Las instalaciones eolicas pequenas evitan los altos costos de extender la 
red electrica a localidades remotas, evitan codes de energi'a y son no contaminantes. De hecho el nivel 
del ruido en el ambiente de la mayoria de turbinas de viento residenciales es de alrededor de 52 a 55 
decibelios. En otras palabras no son tan ruidosas como un refrigerador promedio. 


Datos registrados en estaciones meteorologicas a 10 magi (meters above ground level). Es decir, metros sobre el 
nivel del suelo. 
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Las turbinas de viento convierten la energi'a cinetica del viento en potencia mecanica que mueve a un 
generador para que produzca electricidad limpia. Un sistema basico de turbina aparece en la figura 5-1 
(a) y un sistema eolico aparece en la figura 5-1 (b). 




Figura 5-1 (a) Sistema basico de una turbina eolica y (B) Ejemplo de sistema eolico. 

La mayoria de turbinas fabricadas son maquinas de eje horizontal con dos o tres aspas de fibra de vidrio. 
La cantidad de potencia que una turbina produce depende principalmente del diametro de su rotor. El 
diametro del rotor define el “area de barrido” o la cantidad de viento interceptada por la turbina. En la 
estructura de la turbina se instala el rotor, generador y la cola. La cola mantiene a la turbina frente al 
viento. Debido a que la velocidad se incrementa con la altura, la turbina se monta en una torre. En 
general, entre mas alta sea la torre mas potencia produce el sistema de viento. La formula para 
calcular la potencia de una turbina de viento es : 

P = y 2 k C A p V 3 

P = potencia de salida en Kilowatts 

C = Maximo coeficiente de potencia, es un valor adimensionai que varia de 0.25 a 0.45 (Su maximo valor 

teorico es de 0.59 de acuerdo a la Ley de Betz), 
p = densidad del aire en lb/ft 3 

A = Area de barrido del rotor en ft 2 , tt D 2 /4 , en donde D es el diametro del rotor en ft. 

V = velocidad del viento en mph. 

K = 0.000133 una constante para que el valor de P sea en kilowatts. 

La torre tambien sube a la turbina sobre la turbulencia del aire que existe cerca de la tierra debido a 
obstrucciones tales como montanas, edificios y arboles. Una regia empirica consiste en instalar una 
turbina de viento en una torre con la parte inferior de las aspas del rotor al menos 9 m sobre cualquier 
obstaculo que este a 90 m de la torre (Vea figura 5-2). Subir un generador de 10 KW de una torre de 18 
m de altura a una torre de 30 m implica un incremento del 10% en el costo del sistema, pero se logra 
producir 29% de mas potencia. 
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Figura 5-2 Efecto de los obstaculos sobre la velocidad del viento en un sitio. Observe la ubicacion de las 

aspas de las turbinas para no ser afectadas por la turbulencia. 

Las velocidades del viento se incrementan al subir las alturas de las torres sobre el nivel de la tierra en 
los sitios de viento. Este fenomeno, conocido en la industria del viento como “wind shear”, presenta una 
oportunidad para mejorar la parte economica de una inversion eolica al disponer de turbinas de viento en 
torres mas altas. Por anos, los disenadores de instalaciones de viento han usado la regia empirica de que 
la velocidad del viento en un sitio se incrementa con la ley de potencia de 1/7. La ecuacion matematica 
para esta regia empirica es v 1 = v 0 * (z/z 0 ) a , V 0 es la velocidad del viento a la altura del anemometro, Vi 
es la velocidad del viento a la altura de la torre, Z ^ es la altura de la torre, Z 0 es la altura del anemometro, 
a es el coeficiente de rugosidad con valores tipicos de 1/7, 1/9 y 1/5. 


Ranqos de turbinas eolicas. 

Las pequenas turbinas estan en el rango de 20 W a 100 KW. Las micro turbinas (10 a 500 watts) se usan 
en una variedad de aplicaciones tales como carga de baterias de vehiculos de recreacion y lanchas. Las 
turbinas de 1 a 10 KW se usan en aplicaciones para bombeo de agua. Las turbinas para uso residencial 
varian en tamano de 400 W a 100 KW dependiendo de la electricidad que se quiera generar. Una turbina 
de viento de 1.5 KW llenaria las necesidades de un hogar que requiera de 340 KWh por mes en un lugar 
como en El Salvador con una velocidad media anual de 6.25 m/s. 

El fabricante provee informacion sobre la maxima velocidad a la que la turbina opera adecuadamente. La 
mayoria de turbinas tiene sistemas automaticos de control de sobre velocidad que hacen que el rotor no 
gire fuera de control en situaciones de viento muy altas. 

La potencia en el viento en W/m 2 es el parametro fundamental que determina el exito en un proyecto de 
viento. La produccion de energia del viento a partir de una turbina cambia en proporcion al cubo 
de la velocidad del viento. 


^Como determinar si el potencial eolico de un sitio es adecuado? 


En este documento se sugiere un procedimiento estandar para determinar la capacidad de un sitio para 
ser explotado por un proyecto eolico. El procedimiento tiene 5 pasos y se basa en el embudo mostrado 
en la figura 5-3. 
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Mediciones locales 


Analisis de datos de 
viento (exposicion del 
sitio, rosas de viento, 
etc. 



Embudo para 
estimar el 
potencial de 
viento 




Estimaciones de la velocidad del viento 
en el sitio de interes 


Figura 5-3 Procedimiento para seleccion de un sitio con potencial eolico. 


Los pasos son los siguientes: 

Paso 1. Usar mapas de recursos de viento (mapas eolicos) y tablas de clases de viento. 

La clase de potencia de viento se define como el rango de densidades de potencia del viento (en watts 
por metro cuadrado de area barrida por el rotor, es decir, un area perpendicular al flujo de viento) a una 
altura dado sobre el nivel del suelo. NREL produce estos mapas utilizando un sistema computarizado de 
mapeo analitico y emplrico que combina con modelos matematicos una variedad de conjuntos de datos 
de clima. El mapa eolico de El Salvador desarrollado por NREL en el 2004 aparece en la figura 5-4 y en 
la tabla 5-2 aparecen las clases de viento de NREL 


Tabla 5-2 Definiciones estandar de clases de viento (utilizado en EE.UU. por NREL) 


Clase 

Potencial 

del 

Recurso 

Velocidad del 
viento m/s 
30 m de altura 

Densidad de 
potencia del viento 
W/m 2 

30 m de altura 

Velocidad del 
viento m/s 
50 m de altura 

Densidad de 
potencia del viento 
W/m 2 

50 m de altura 

1 

Pobre 

0.0 -5.1 

0-160 

0-5.6 

0-200 

2 

Marginal 

5.1 -5.9 

1 60 - 240 

5.6 -6.4 

200 - 300 

3 

Considerable 

5.9 -6.5 

240 - 320 

6.4 -7.0 

300 - 400 

4 

Bueno 

6.5 -7.0 

320 - 400 

7.0 -7.5 

400 - 500 

5 

- 

7.0 -7.4 

400 - 480 

7.5 -8.0 

500 - 600 

6 

- 

7.4 -8.2 

480 - 640 

8.0 -8.8 

600 - 800 

7 

- 

8.2-11.0 

640-1600 

8.8-11.9 

800 - 2000 


Nota : La densidad de potencia se expresa en watts por metro al cuadrado (W/m z ) de area barrida 
perpendicular al flujo del viento. 
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Figura 5-4 Mapa eolico de El Salvador elaborado por NREL. 


Paso 2. Indirectamente cuantificar el recurso de viento a partir de informacion de viento de estaciones 
meteorologicas cercanas. 

Datos deseables de mediciones de viento cercanas al sitio de interes son: 

■ Elevacion del sitio. 

■ Promedio mensual de la velocidad del viento. 

■ Rosa de viento (dato de la frecuencia de la velocidad y direccion del viento, que muestran 
direcciones en las que el viento es mas fuerte y mas debil). Las rosas de viento son una especie 
de huella meteorologica. Un vistazo a la rosa de vientos es extremadamente util para situar 
aerogeneradores (turbinas eolicas). 

■ Exposicion del sitio (terreno local- montanas, valles, etc.- cobertura vegetativa y estructuras 
construidas por el hombre). Un sitio bien expuesto es aquel que no tiene influencias de arboles o 
edificios. 

■ Rugosidad de la superficie (ver tabla 5-3). 

■ Alturas sobre el nivel del suelo. 

■ Recuperacion de datos (numero de horas de datos validos versus total de horas posibles). 

■ Localizacion del sitio con respecto a su propiedad. 

Una vez validados los datos los puede usar para estimar el recurso viento en su propiedad. 


Tabla 5-3 Tabla de clases y de longitudes de rugosidad. 


Clase de 
rugosidad 

Longitud de 
rugosidad en m 

Tipo de paisaje 

0 

0.0002 

Superficie del agua 

0.5 

0.0024 

Terreno completamente abierto con una superficie lisa, p.ej., pistas 
de hormigon en los aeropuertos, cesped cortado, etc. 

1 

0.03 

Area agricola abierta sin cercados ni setos y con edificios muy 
dispersos. Solo colinas suavemente redondeadas . 

1.5 

0.055 

Terreno agricola con algunas casas y setos resguardantes de 8 
metros de altura con una distancia aproximada de 1250 m. 
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Clase de 
rugosidad 

Longitud de 
rugosidad en m 

Tipo de paisaje 

2 

0.1 

Terreno agricola con algunas casas y setos resguardantes de 8 
metros de altura con una distancia aproximada de 500 m. 

2.5 

0.2 

Terreno agricola con muchas casas, arbustos y plantas, o setos 
resguardantes de 8 metros de altura con una distancia aproximada 
de 250 m. 

3 

0.4 

Pueblos, ciudades pequehas, terreno agricola, con muchos o altos 
setos resguardantes, bosques y terreno accidentado y muy 
desigual 

3.5 

0.8 

Ciudades mas grandes con edificios altos 

4 

1.6 

Ciudades muy grandes con edificios altos y rascacielos 


Nota: Definiciones de acuerdo con el Atlas Eolico Europeo, WAsP. 


Paso 3. Otra manera indirecta de cuantificar el recurso es observar la vegetacion del area. Los arboles, 
especialmente, las coniferas pueden deformarse permanentemente por vientos fuertes. Esta deformacion 
se conoce como “flagging” y ha sido utilizada para estimar la velocidad del viento promedio en un area 
(ver figura 5-5). 


Flagging 


Prevailing wind 



7777777777 


II 

Slight 

flagging 



0 

No 

deformity 


7777777777 



I 

Brushing 
and slight 
flagging 


7777777777 




III 

Moderate 

flagging 


7777/77777 



IV 

Complete 

flagging 


7J77Tm 



VI 

Complete 

throwing 


VII 

Carpeting 



02979310m 


Griggs-Putnam Index of Deformity 

Index I II III IV V VI VII 

Wind mph 7-9 9-11 11-13 13-16 15-18 16-21 22+ 

Speed tn/s 3-4 4-5 5-6 6-7 7-8 8-9 10 


Figura 5-5 Indice de deformidad (“flagging”) de Griggs-Putnam. 



Figura 5-6 Ejemplos de “flagging”, (a) Clase II correspondiente a velocidades de 4 a 5 m/s; (b) Clase VI 
correspondiente a velocidades de 8 a 9 m/s. 
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Paso 4. Monitoreo directo por un sistema de medicion del recurso de viento da una idea clara del recurso 
disponible. 

Las estimaciones mas confiables se obtienen a traves de mediciones en el sitio especffico. Esto puede 
ser un poco costoso y requiere de 6 meses a varios anos para completarse. 

Paso 5. Ademas de medir o estimar la velocidad media anual en el sitio se necesita conocer las 
direcciones de viento que mas prevalecen en su sitio. 

El conjunto ideal de datos para un proyecto de viento consistirfa de una decada de mediciones de 
velocidad y direccion del viento junto con la temperatura y presion tomada sobre, abajo y a la altura 
exacta (sobre la superficie de la tierra) de la torre de la turbina de viento seleccionada. 


Distribucion de Weibull v Rosa de los vientos. 


Para la industria eolica es muy importante ser capaz de describir la variacion de las velocidades del 
viento. Los proyectistas de turbinas necesitan la informacion para optimizar el diseno de sus 
aerogeneradores, asi como para minimizar los costos de generacion. Los inversores necesitan la 
informacion para estimar sus ingresos por produccion de electricidad. La variacion del viento en un 
emplazamiento tipico (estacion meteorologica) suele describirse utilizando la llamada Distribucion de 
Weibull, como la mostrada en la figura 5-7 para el caso de la estacion A35 del Cerro Verde. Esta 
estacion en particular tiene una velocidad media del viento de 3.91 metros (parametro U) por segundo, y 
la forma de la curva esta determinada por un parametro de forma de 1.32 (parametro k). El grafico 
muestra una distribucion de probabilidad. El area bajo la curva siempre vale exactamente 1, ya que la 
probabilidad de que el viento sople a cualquiera de las velocidades, incluyendo el cero, debe ser del 100 
%. 


Como se podra observar, la distribucion de las velocidades del viento es sesgada, es decir, no es 
simetrica. A veces tendra velocidades de viento muy altas, pero son muy raras (en este caso de 6.17 
m/s). Por otro lado, las velocidades del viento de 1.5 m/s son las mas comunes. Los 1.5 metros por 
segundo es el llamado valor modal de la distribucion. Si multiplicamos cada diminuto intervalo de la 
velocidad del viento por la probabilidad de tener esa velocidad particular, y los sumamos todos, 
obtenemos la velocidad del viento media. 

La distribucion estadistica de las velocidades del viento varia de un lugar a otro del globo, dependiendo 
de las condiciones climaticas locales, del paisaje y de su superficie. Por lo tanto, la Distribucion de 
Weibull puede variar tanto en la forma como en el valor medio. Si el parametro de forma es exactamente 
2 la distribucion es conocida como distribucion de Rayleigh. En el caso de El Salvador la estacion mas 
cercana a una distribucion de rayleigh es Z-4 la Galera en Morazan que tiene un factor de forma de 1 .73. 

La rosa de los vientos muestra informacion sobre las distribuciones de velocidades del viento y la 
frecuencia de variacion de las direcciones del viento. La rosa se divide en 12 sectores cada uno de 30° 
del horizonte. Un vistazo a la rosa de los vientos es extremadamente util para situar un aerogenerador. Si 
una gran parte de la energia del viento viene de una direccion en particular la turbina eolica se ubicara en 
esa direccion, asegurando que en esa direccion exista la menor cantidad de obstaculos y que el terreno 
sea lo mas liso posible. 
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'A-35' Observed Wind Climate 


m 


H n a™ \i I ii klfi 


swera 


Produced on IS/ 06/ 200 5 st 06:54:09 a nn by licenced user: Unlicenced using WAsP Version : 
8 01,0057 


Site description: 'Cerro Verde'; position: 13.B3°N -S9,62 n E; anemometer height: 10 00 magi 


- 

Unit 

Measured 

Weibull Fit 

Discrepancy 

Mean wind speed 

m/s 

unknown 

3.91 

u n kn own 

Mean power density 

W/m z 

unknown 

123 

u n kn own 



- 

O 

30 

60 

90 

120 

150 

ISO 

210 

240 

270 

300 

330 
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A 

6.S 

5.2 

0 5 

2.0 

0.5 

1.8 

1.7 

1.8 

0.5 

1.1 

0.5 

2.7 

4.2 

k 

1.49 

1.63 

10.34 

2.00 

10.34 

2.16 

2.02 

1.42 

10.34 

1.24 

10.34 

0.77 

1.32 

y 

6.17 

4.62 

0.49 

1.74 

0.49 

1.58 

1=51 

1.67 

0.49 

1=02 

0.49 

3.14 

3.91 

p 

397 

146 

0 

6 

0 

4 

4 

9 

0 

2 

0 

241 

123 

Freq 

2.9 

72.4 

0.2 

3.4 

0.2 

4.0 

2.8 

10.1 

0.2 

1.8 

0.2 

1.8 

100 


Figura 5-7 Huella meteorologica de la estacion A-35 Cerro Verde, Santa Ana. Muestra la rosa de los 

vientos a la izquierda y la distribucion de Weibull a la derecha. 


5.2 Clima eojjco de las estaciones meteoroloqicas del pais. 


La base de datos del clima eolico del proyecto SWERA-UCA posee un total de 20890 horas de registros 
de velocidad y direccion del viento en 18 estaciones meteorologicas, ver figura 5-8. Se encuentra 
compuesta de datos digitalizados a partir de los cuadros mensuales de anemometros de la division de 
Meteorologia de SNET, y corresponden a los registrados durante los anos de 1980-1986. Cada registro 
corresponde a la medicion horaria de la velocidad y direccion del viento durante un dia de medicion mas 
el valor de direccion y hora en la que ocurrio la maxima velocidad instantanea del dia. 

Los datos se encuentran almacenados en una base de datos de Microsoft Access y basicamente se 
pueden describir en dos de sus tablas: medidas y detalles de las medidas, como se muestra en la figura 
5-9. La primera tabla registra la fecha de medicion y por medio de la Have primaria “banda” se relaciona 
con la segunda tabla. Por ejemplo el registro de viento de la estacion L4 del dia 18 de marzo de 1983 le 
corresponde una Nave primaria 17058 y buscando este numero en la segunda tabla se encuentran tres 
registros, el primero de la direccion del viento en decenas de grados, el segundo de velocidad en km/h 
para cada una de las horas del dia y el tercero el valor de la maxima velocidad instantanea del viento 
(rafaga) detectada en el dia. Para este mismo dia (18 de marzo) le corresponde a las 10 de la manana 
(columna hi 0) una direccion de viento predominante de 28 decenas de grado, es decir 280 grados de 
acimut y una velocidad de 40 km/h. La maxima velocidad instantanea registrada ese dia fue de 6.9 m/s 
a las cuatro de la tarde (columna hi 6) y corresponde a una direccion de viento predominante de 22 
decenas de grados (columna rumbo_r) que equivale a un acimut de 220 grados. 
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Figura 5-8 Estaciones meteorologicas que poseen registros horarios de viento. 



Figura 5-9 Estructura de las tablas principales de la base de datos. 


Para la determinacion de la velocidad y la potencia media mensual, anual y por temporada (estaciones 
climaticas: lluviosa, seca y sus transiciones) se utilizo el calculo de la densidad de probabilidad por medio 
de la funcion de Weibull utilizando el programa WAsP. Los datos para el calculo en WAsP fueron 
obtenidos por medio de consultas SQL a la base de datos de viento de tal manera que se dispusiera de 
pares ordenados de velocidades y direcciones para cada hora durante periodos de tiempo determinados. 
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Llagandose a obtener el clima eolico en cada una de las estaciones meteorologicas para cada mes, ano y 
temporada de cada una de las 18 estaciones meteorologicas. Los valores calculados corresponden a: 
factores de forma y escala de Weibull, velocidad media y potencia media (utilizando la funcion) y 
distribucion de velocidades por sector de medicion en porcentaje en base al miliar. En la figura 5-10 se 
muestra el clima eolico anual de la estacion S-10 en la cual la velocidad media fue de 2.44 m/s y la 
potencia 22 W/m 2 , adicionalmente se espera una potencia de 30 W/m 2 del viento procedente del sector 1 
es decir del Norte. 


'S-10' Observed Wind Climate 

Produced on 16/06/2005 at OS: 19: 14 a.m. by licenced user: Ismael Antonio Sand"iez r Universidad 
Centroamericana, San Salvador using WAsP Version: B .01.00 57 


Site description: 'Aeropuerto Ilopango 1 ; position: 13.70°N -89.12 n E; anemometer height: 10.00 m 
a gJ, 


- 
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P 
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17 

22 
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Figura 5-10 Datos principales del clima eolico de la estacion meteorologica, calculado en WAsP. 


5.2.1 Climatoloqia regional 


A partir del clima eolico e informacion de la topografia de aquellos lugares que presentaran los mayores 
valores de potencias medias se determino la climatologia regional a 50 m.a.g.l. Como se muestra en la 
figura 5-1 1 los sitios de mayor potencial eolico fueron: Cerro Verde, Planes de Montecristo y La Galera. 
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Figura 5-1 1 Potencia media registrada, la primera columna corresponde a la temporada seca, la segunda 

a la seca-lluviosa, la tercera a la lluviosa y la cuarta a la lluviosa seca . 

La extension territorial alrededor de cada uno de los lugares evaluados fue aproximadamente de cuadros 
de 10 km por lado centrados en la ubicacion de las estaciones, a excepcion de la estacion que registro el 
mayor potencial eolico (Cerro Verde) que se tomo un sector mayor con el objeto de abarcar al lago de 
Coatepeque. 

La rugosidad del lugar fue calculada a partir del mapa de uso de suelo de El Salvador de 1996, 
publicado en el Sistema de Informacion Ambiental de El Salvador; al reclasificar los diferentes usos de 
suelo en la zona. El mapa de uso de suelo fue creado al georeferenciar y clasificar las bandas 3, 4 y 5 de 
imagenes Landsat - TM tomadas entre 1993 y 1994; trabajo realizado por el Centro agronomico tropical 
de investigacion y ensenanza CATIE y la direccion de economia agropecuaria DGEA del Ministerio de 
Agricultura y Ganaderia MAG. Los valores utilizados para la rugosidades son presentados en la tabla5-4 
y corresponden a valores recomendados por Humberto Rodriguez 2 en el proyecto de Implementacion de 
medidas para el desarrollo de un mercado nacional sostenible para servicios de energia en el area rural. 

Tabla 5-4 Asignacion de rugosidades en base a los tipo de uso de suelo usado para el modelaje eolico 

en WAsP. 


Uso de suelo 

Rugosidad 

Areas urbanas 

0.32 

Bosques naturales 

0.24 

Cafe 

0.24 

Cana de azucar 

0.18 

Cuerpo de agua 

0.1 

Lava 

0.14 

Pastos y granos 
basicos 

0.1 


Consultor colombiano del MINEC en los temas solar y eolico. 
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Figura 5-12 Zonas seleccionadas para evaluar la climatologfa regional 


La orografia del lugar fue obtenida por medio de la digitalizacion de las curvas de nivel a cada 20 m de 
cartas topograficas a escala 1:50,000. Las cartas fueron escaneadas y georeferenciadas para su 
posterior digitalizacion en el laboratorio de sistemas de informacion geografica de la UCA. El 
procedimiento seguido para la creacion de las capas tematicas consistio en el escaneo de mapas a 
escala 1:50,000 publicados por el Centro Nacional de Registros y su posterior georreferenciacion. 
Posteriormente se procedio a realizar una digitalizacion de las curvas de nivel y de los puntos de 
elevacion conocidos. Este trabajo fue realizado en el sistema de informacion geografico ArcGIS de ESRI. 

La rugosidad y la orografia obtenida en formatos de sistemas de informacion geograficas fueron 
convertidas a archivos MAP de WAsP por medio de la modificacion de un Script programado en Avenue 
por Javier Gonzalez y publicado en la web de los fabricantes del software de sistemas de informacion 
geografica, ArcView. 

Los resultados obtenidos del modelaje en WAsP consistieron en informacion de la funcion de Weibull 
para cada localidad del area de estudio y para cada sector de direccion de viento. Se exportaron a un 
archivo ASCII de extension WRG (ver la tabla 5-5 para el contenido del archivo generado) para su 
posterior incorporacion a los sistemas de informacion geografica y comparacion con el atlas eolico de 
Centro America. 


Tabla 5-5 Descripcion del archivo WRG de intercambio. 


Columns 

Width 

Description 

01-10 

10 

GridPoint 

11-20 

10 

X-coordinate (easting) of the site [m] 

21-30 

10 

Y-coordinate (northing) of the site [m] 

31-38 

8 

Z-coordinate (elevation) of the site [m] 

39-43 

5 

Height above ground level [m a.g.l.] 

44-48 

5 

Weibull A-parameter for the total distribution [ms 1 ] 

49-54 

6 

Weibull /c-parameter for the total distribution 

55-69 

15 

Power density [Wm 7 ] or power production [Why 1 ] 
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Columns 

Width 

Description 

70-72 

3 

Number of sectors 

73-76 

4 

Frequency of occurrence for sector #1 [%■ 1 0] 

77-80 

4 

Weibull A-parameter for sector #1 [ms 1 -10] 

81-85 

5 

Weibull /c-parameter for sector #1 [-100] 

86-98 

13 

As columns 73-85, but for sector #2 

■ ■ ■ 

■ ■ ■ 

■ ■ ■ 

216-228 

13 

As columns 73-85, but for sector #12 


5.3 Climatoloqla Regional en el Cerro Verde (A-35) 


Descripcion general 

La region de estudio en Cerro Verde consiste en un cuadro de 16 km de lado que se encuentra, 
aproximadamente, delimitado por las faldas del volcan de Santa Ana e Izalco al norte y oeste, por la 
carretera de Sonsonate al sur y por el li'mite de la cuenca de Coatepeque al este. En el cerro verde se 
encuentra una estacion meteorologica con 978 registros de viento diarios y utilizando la funcion de 
Weibull se ha determinado valores de velocidades y potencias medias anuales de 3.91 m/s y 123 W/m 2 
al ano con una direccion de viento predominante NNE, en la figura 5-13 se presenta, por medio de la rosa 
de vientos, la direccion de viento predominante y la distribucion de frecuencias de las velocidades y los 
factores de la funcion de Weibull. En las tablas 5-6 y 5-7 se muestran los factores de Weibull, velocidad y 
potencia media para otros periodos de tiempo. En los periodos de tiempo evaluados la direccion 
predominate del viento es siempre NNE. 




Figura 5-13 Rosa de los vientos y distribucion de velocidades para la estacion A-35 (Cerro Verde). 
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Tabla 5-6 Factores de Weibull, Velocidad y Potencia media para cada uno de los meses del ano. 


Parametros 

Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

A 

5 


3.3 

3.8 

2.5 

2.8 

5.6 

4.4 

4 

4.7 

5.7 

6 

K 

1.46 

1.26 

1.21 

1.24 

1.22 

1.1 

2.04 

1.61 

1.83 

1.4 

1.6 

1.59 

U 

4.55 

4.76 

3.09 

3.55 

2.35 

2.67 

5 

3.95 

3.55 

4.28 

5.08 

5.38 

P 

163 

241 

71 

103 

31 

56 

143 

93 

58 

145 

200 

238 


Parametros: A, factor de escala; K, factor de forma; U, velocidad (m/s); P, potencia (W/m 2 ) 

Tabla 5-7 Factores de Weibull, Velocidad y Potencia media para cada temporada del ano. 


Parametros 

Temporada 

Anual 

1 

2 

3 

4 

A 

4.7 

2.6 

4.2 

5 

4.2 

k 

1.29 

1.31 

1.54 

1.36 

1.32 

U 

4.31 

2.41 

3.8 

4.57 

3.91 

P 

170 

29 

88 

186 

123 


Parametros: A, factor de escala; K, factor de forma; U, velocidad (m/s); P, potencia (W/m 2 ) 
Temporada: s, seca; s-ll, seca lluviosa; II, lluviosa y ll-s, lluviosa seca. 


Ruaosidad 

La rugosidad de la region vari'a de 0.1 en el lago hasta 0.32 en zonas urbanas, siendo el valor 
predominante en el 77% de la region la rugosidad de 0.24, asignada a los bosques naturales y los 
cafetales. La estacion de medicion se encuentra ubicada en el cerro verde en una zona de vegetacion 
espesa con valores de rugosidad de 0.24, as! como se muestra en la figura 5-14. 



Figura 5-14 Mapa de uso de suelo reclasificado con valores de rugosidad en los alrededores de la 

estacion Cerro Verde (A-35). 
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Oroqrafia 


El relieve en la zona es de tipo montanoso con alturas maximas de 2370 y rmnimas de 560 msnm en 
distancias cortas. Existen tres puntos de mayor elevacion en la zona de los cuales dos de ellos son 
volcanes y el otro es un cerro. Los volcanes dentro de la zona son el Lamatepec o volcan de Santa Ana y 
el volcan de Izalco. El cerro es conocido como Cerro Verde y es en el que se encuentra ubicada la 
estacion meteorologica a una altura de 2030 msnm Al este de la zona montanosa se ubica el lago de 
coatepeque a una altura de 750 msnm, rodeado el lago en casi toda su totalidad por una cordillera con 
amplia vegetacion, verfigura 5-15. 



Figura 5-15. Orografia de la region “Cerro Verde”, (a) Curvas de nivel de la zona; (b) Modelo de elevacion 
digital. 
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Estimacion del potencial eolico de la zona Cerro Verde 


La velocidad media de la region a 50 m.a.g.l. es mostrada en la figura 5-16, identificandose cuatro lugares 
con las magnitudes mayores de velocidades: Cerro Verde y Cerro Mala Cara, Volcanes de Izalco y de 
Santa Ana. Las menores velocidades de viento se presentan en zonas bajas o en aquellas que el efecto 
topografico les ha afectado grandemente. La figura 5-17 muestra el potencial eolico a 50 m.a.g.l. de la 
zona “Cerro Verde” obtenido por medio del modelo realizado en WAsP y el atlas eolico de Centro 
America. La diferencia en las magnitudes de potencia se puede atribuir a que el modelaje realizado en 
WAsP utiliza datos a 10 m.a.g.l. de estaciones de larga duracion, por lo que es posible obtener en dicha 
zona una mayor potencia energetica que la calculada. Es importante senalar la existencia de sitios 
energeticos en el Cerro Mala Cara y la zona del cerro Mascaron que no quedaban definidos en el Atlas 
eolico deNREL. 



Figura 5-16 Clima regional de la zona “Cerro Verde” (Velocidades a 50 m.a.g.l.). 
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Cerro Mala Cara 




^Mnirrmidt- 

m • 

Hnawii\il j 




ri 


swera 



2 . 3 - 57.2 
57.3 - 112.1 
1 12.2 - 157.0 
157.1 - 221.9 
222.0 - 276.8 


Figura 5-17 Mapas de potencial eolico a 50 m.a.g.l. en W/rrr obtenido por medio del modelaje de la 

region “Cerro Verde”. 


5.4 Climatoloqla Regional en Planes de Montecristo (A-31). 


Descripcion general 


La region de estudio consiste en un cuadro de 10 km de lado que se encuentra, aproximadamente, 
delimitado por la frontera de Honduras al norte, por la quebrada la Chicotera al este y por la carretera que 
conduce de Metapan a la Hacienda El Pinar al sur y por los cerros Copacho, Carrapinga y Anona. Al 
centra de la region, en la Loma Las Toreras se encuentra una estacion meteorologica (Planes de 
Montecristo) con 1044 registros de viento diarios y utilizando la funcion de Weibull se han determinado 
valores de velocidades y potencias medias anuales de 2.65 m/s y 45 W/m 2 al ano con una direccion de 
viento predominante NNE todo el ano y en algunos meses SSW y E, en la figura 5-18 se presenta por 
medio de la rasa de vientos la direccion de viento predominante y la distribucion de frecuencias de las 
velocidades y los factores de la funcion de Weibull. En las tablas 5-8 y 5-9 se muestran los factores de 
Weibull, velocidad y potencia media para otros periodos de tiempo. 
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Figura 5-18 Rosa de los vientos y distribution de velocidades para la estacion A-31 . 


Tabla 5-8 Factores de Weibull, Velocidad y Potencia media para cada uno de los meses del afio. 



Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

A 

3 

4.2 

2.6 

3 

1.8 

1.7 

3.2 

2.9 

1.9 

3.3 

4 

4 

k 

1.23 

1.44 

1.21 

1.15 

1.41 

1.03 

1.41 

1.41 

1.39 

1.45 

1.53 

1.44 

U 

2.8 

3.83 

2.4 

2.83 

1.64 

1.71 

2.92 

2.65 

1.77 

3 

3.57 

3.61 

P 

51 

100 

33 

60 

8 

17 

46 

34 

10 

47 

74 

83 


Parametros: A, factor de escala; K, factor de forma; U, velocidad (m/s); P, potencia (W/m 2 ) 


Tabla 5-9 Factores de Weibull, Velocidad y Potencia media para cada temporada del afio. 



Temporada 

Anual 

1 

2 

3 

4 

A 

3.4 

1.9 

2.3 

4.4 

2.8 

k 

1.28 

1.33 

1.2 

1.91 

1.21 

U 

3.12 

1.78 

2.17 

3.93 

2.65 

P 

66 

11 

25 

74 

45 


Parametros: A, factor de escala; K, factor de forma; U, velocidad (m/s); P, potencia (W/m 2 ) 
Temporada: s, seca; s-ll, seca lluviosa; II, lluviosa y ll-s, lluviosa seca. 


Ruqosidad 


La rugosidad de la region es de dos tipos: 0.1 en pastos y granos basicos que cubre el 25% de la region y 
de 0.24 a los bosques naturales y a los cafetales cubriendo el resto del territorio. La estacion de medicion 
se encuentra ubicada en Loma las Toreras, zona de vegetacion espesa con valores de rugosidad de 
0.24, asi como se muestra en la figura 5-19. 
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Figura 5-19 Mapa de uso de suelo reclasificado con valores de rugosidad. 


Orografia 

El relieve en la zona es de tipo montanoso con alturas maximas de 2400 y minimas de 800 msnm en 
distancias cortas, como se muestra en la figura 5-20. La estacion meteorologica se encuentra a una 
altura de 1831 msnm. La region se encuentra en una zona montanosa con grandes pendientes en cortas 
distancias, con abundante vegetacion de tal manera que la parte oeste es una zona nacional protegida 
llamada “Montecristo”. 



Figura 5-20 Orografia de la region Planes de Montecristo (Metapan), (a) Curvas de nivel; (b) Modelo de 
elevacion digital. 
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Estimacion del potencial eolico de la zona Planes de Montecristo. 


La velocidad media de la region a 50 m.a.g.l. es mostrada en la figura 5-21, identificandose cinco lugares 
con las magnitudes mayores de velocidades: la zona de Miramundo y El Cerro Los Cantaros (ubicados 
en la zona protegida), Cerro Redondo, Loma El Malcotal y Cerro Guamilar. Las menores velocidades de 
viento se presentan en zonas bajas o en aquellas que el efecto topografico les ha afectado grandemente. 
La figura 5-22 muestra el potencial eolico a 50 m.a.g.l. de la zona, obtenido por medio del modelo 
realizado en WAsP y el atlas eolico de NREL para Centroamerica. Existen zonas al sur de la region que 
el Atlas de NREL las establece con alto potencial energetico que no aparecen en el modelo realizado, en 
las cuales se requiere de ubicar una estacion de corta duracion en la zona para evaluar su potencial 
energetico. La diferencia en las magnitudes de potencia se puede atribuir a que el modelaje realizado en 
WAsP utiliza datos a 10 m.a.g.l. de estaciones de larga duracion, por lo que es posible obtener en dicha 
zona una mayor potencia energetica que la calculada. 



Figura 5-21 Clima regional de la zona Planes de Montecristo (Velocidades a 50 m.a.g.l.) 



Figura 5-22 Mapas de potencial eolico a 50 m.a.g.l. en W/m 2 obtenido por medio del modelaje de la 

region Planes de Montecristo. 
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5.5 Climatoloqia Regional en La Galera (Z-4). 


Descripcion general 

La region de estudio consiste en un cuadro de 10 km de lado que se encuentra, aproximadamente, 
delimitado por la frontera de Honduras al norte, por el Cerro Bailaderos y Las Cumbres del Cerro Verde al 
este y por los Cerros Las Crucitas y El Colorado al sur y por los cerros El Infernal y La Ardilla. Al centro 
de la region, en Las Galeras se encuentra una estacion meteorologica con 618 registros de viento diarios 
y utilizando la funcion de Weibull se ha determinado valores de velocidades y potencias medias anuales 
de 2.88 m/s y 33 W/m 2 al ano con una direccion de viento predominante N y NNE todo el ano, en algunos 
meses las direcciones son variadas con velocidades pequenas, en la figura 5-23 se presenta por medio 
de la rosa de vientos la direccion de viento predominante y la distribucion de frecuencias de las 
velocidades y los factores de la funcion de Weibull. En las tablas 5-10 y 5-1 1 se muestran los factores de 
Weibull, velocidad y potencia media para otros periodos de tiempo. 
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Figura 5-23 Rosa de los vientos y distribucion de velocidades para la Z-4. 

Tabla 5-10 Factores de Weibull, Velocidad y Potencia media para cada uno de los meses del ano. 



Ene 

Feb 

Mar 

Abr 

May 

Jun 

Jul 

Ago 

Sep 

Oct 

Nov 

Die 

A 

4 

5.1 

2.7 

3 

2.1 

2.8 

3.5 

3.6 

2.6 

3.7 

4.3 

4.9 

k 

1.65 

2.16 

1.4 

1.74 

1.54 

1.83 

2.41 

2.24 

2.16 

1.89 

1.99 

2.44 

U 

3.54 

4.54 

2.5 

2.66 

1.85 

2.47 

3.14 

3.21 

2.26 

3.28 

3.79 

4.34 

P 

65 

102 

29 

26 

10 

20 

31 

35 

13 

44 

64 

80 


Parametros: A, factor de escala; K, factor de forma; U, velocidad (m/s); P, potencia (W/m 2 ) 


Tabla 5-1 1 Factores de Weibull, Velocidad y Potencia media para cada temporada del ano. 



Temporada 

Anual 

1 

2 

3 

4 

A 

3.8 

2.3 

3 

4.3 

3.2 

k 

1.72 

1.68 

1.97 

1.95 

1.73 

U 

3.35 

2.06 

2.66 

3.79 

2.88 

P 

52 

13 

22 

65 

33 
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Parametros: A, factor de escala; K, factor de forma; U, velocidad (m/s); P, potencia (W/m 2 ) 
Temporada: s, seca; s-ll, seca lluviosa; II, lluviosa y ll-s, lluviosa seca. 




Ruqosidad 


Segun el mapa de uso de suelos utilizados el tipo de suelo de la region es Cafetal o Bosques Naturales 
por lo que se utilizo un factor de rugosidad de 0.24 en toda la zona. 

Oroqrafi'a 


El relieve en la zona es de tipo montanoso con alturas maximas de 2020 y minimas de 1120 msnm en 
distancias cortas, como se muestra en la figura 5-24. La estacion meteorologica se encuentra a una altura 
de 1900 msnm. La region se encuentra en una zona montanosa con grandes pendientes en cortas 
distancias, con abundante vegetacion. 



Figura 5-24 Orografia de la region “La Galera”. (a) Curvas de nivel; (b) Modelo de elevacion digital. 


Estimacion del clima regional de la zona La Galera 


La velocidad media de la region a 50 m.a.g.l. es mostrada en la figura 5-25, identificandose tres lugares 
con las magnitudes mayores de velocidades: La Loma de Enmedio y Costado Nor-Este de Bailaderos y el 
Volcan de Agua. Las menores velocidades de viento se presentan en zonas bajas o en aquellas que el 
efecto topograficos les ha afectado grandemente. La figura 5-26 muestra el potencial eolico a 50 m.a.g.l. 
de la zona, obtenido por medio del modelo realizado en WasP y el atlas eolico de NREL para 
Centroamerica. La ausencia de datos del Atlas se debe a los h'mites politicos de El Salvador, pero se 
observa similitud en la distribution de las potencias y no a la magnitud de las mismas. La diferencia en 
las magnitudes de potencia se puede atribuir a que el modelaje realizado en WAsP utiliza datos a 10 
m.a.g.l. de estaciones de larga duration, por lo que es posible obtener en dicha zona una mayor potencia 
energetica que la calculada. 
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Figura 5-25 Clima regional de la zona “La Galera” (Velocidades a 50 m.a.g.l.) 



Figura 5-26 Mapas de potencial eolico a 50 m.a.g.l. en W/m 2 obtenido por medio del modelaje de la 

region La Galera. 
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5.6 Tablas de medicion de potencial eolico en las estaciones 

meteoroloqicas del pais. 


El proyecto SWERA ejecutado por la UCA considero datos de viento horarios de 18 estaciones las 
cuales se presentan en la tabla 5-12. El total de datos para modelaje en WAsP utilizo 20,890 horas de 
mediciones de todas las estaciones. Sus resultados en W/m 2 y KWh/m 2 se presentan en la tabla 5-13. 


Tabla 5-12 Estaciones con mediciones de viento consideradas en el proyecto SWERA. 


Estacion 

Nombre 

Longitud 

Latitud 

Elevacion 

A12 

Santa Ana El Palmar 

-89.5700 

13.9766 

725 

A31 

Planes de Montecristo 

-89.3600 

14.3983 

1851 

A35 

Cerro Verde 

-89.6233 

13.8266 

2030 

BIO 

Cerron Grande 

-88.7850 

13.9383 

600 

G3 

Nueva Concepcion 

-89.2900 

14.1250 

320 

H13 

Apaneca Santa Leticia 

-89.7816 

13.8616 

1300 

H8 

Ahuachapan 

-89.8600 

13.9433 

725 

L18 

El Boqueron 

-89.2816 

13.7350 

1800 

L4 

San Andres 

-89.4067 

13.8083 

460 

M23 

Cerro Cacahuatique 

-88.2317 

13.7633 

1975 

M6 

El Papalon 

-88.1267 

13.4400 

1960 

N15 

La Union 

-87.8816 

13.3316 

95 

S10 

Aeropuertpo-llopango 

-89.1183 

13.6983 

615 

S27 

Estacion Matriz 

-89.1433 

13.6883 

650 

U11 

Beneficio La Carrera 

-88.5250 

13.3300 

75 

VI 3 

Apastepeque 

-88.7466 

13.7000 

570 

V6 

Santa Cruz Porrillo 

-88.8033 

13.4400 

30 

Z4 

La Galera 

-88.0866 

14.0430 

1900 


En la tabla 5-14 se presenta una comparacion entre los datos de viento obtenidos por la UPES en 1987 y 
los obtenidos por el proyecto SWERA de la UCA utilizando WAsP. 

Tabla 5-4 Areas que presentan mayor potencial eolico de acuerdo a estudio UPES y al proyecto SWERA- 

UCA (calculados con WAsP). Datos a 10 magi. 


Codigo 

Nombre y 
depto. 

Velocidad 

media 

anual 

(m/s) 

UPES 

Velocidad 

media 

anual 

(m/s) 

SWERA 

Densidad 

de 

potencia 

media 

anual 

(W/m 2 ) 

UPES 

Densidad 

de 

potencia 

media 

anual 

(W/m 2 ) 

SWERA 

Densidad 
de energia 
media 
anual 
(KWh/m 2 ) 
UPES 

Densidad 
de energia 
media 
anual 
(KWh/m 2 ) 
SWERA 

A-35 

Cerro Verde 
(Santa Ana) 

4.39 

3.91 

116 

123 

1040 

1077 

A-31 

Planes de 
Montecristo 
(Santa Ana) 

3.14 

2.65 

70 

45 

629 

394 

Z-4 

La Galera 
(Morazan) 

2.86 

2.88 

57 

33 

519 

289 
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Tabla 5-3 Datos calculados por WAsP para cada mes del ano, se presentan datos en W/m 2 y KWh/m 2 (anual). 



lenuy ojpeiuojd 

105 

394 

1077 

70 

70 

CD 

CD 

70 

193 

53 

149 

53 

140 

193 

44 

70 

123 

CO 

289 

CM 

45 

123 

00 

00 

03 

00 

22 

CD 


CD 

CD 

22 

LO 

00 


1^ 

33 

0 jqLU 0 P!Q 

158 

727 

2085 

CO 

140 

473 

114 

420 

CO 

237 

CO 

184 

324 

53 

CO 


00 

00 

701 

00 

83 

238 


CD 

54 

00 

00 


27 

1^ 

CM 

37 

CD 



o 

o 

00 

0 jqiu 0 jAON 

00 

00 

00 

CD 

1752 

35 

CO 

315 

70 

394 

44 

123 

26 

79 

219 

35 

53 


44 

5 

LO 

o 

74 

200 



CD 

00 

00 

45 

LO 


CO 

03 

25 


CD 


LO 

co 

©jqnpo 

70 

412 

1270 


26 

403 

35 

184 

00 

123 

03 

44 

123 

26 

44 

53 

26 

385 

00 

47 

145 


00 

46 


CM 

CM 


- 

LO 


CO 

LO 

CD 

00 

PP 

0 jqiu 0 i;d 0 S 

44 

00 

00 

508 

35 

03 

44 

CD 

CM 

00 

00 

26 

00 

00 

00 

35 

79 

26 

44 

00 

00 

00 

114 

LO 

o 

58 


- 

LO 

00 

o 

CO 

o 

CM 


03 

CO 

LO 

O 

CM 

00 

ojsoBy 

53 

298 

815 

44 

26 

79 

CD 

CM 

70 

26 

140 

26 

44 

00 

00 

26 

53 

96 

35 

307 

CD 

34 

93 

LO 

00 

03 

00 

00 

CO 

CD 

CO 

LO 

O 

CO 

CD 

- 

^f- 

35 

oi|np 

53 

403 

1253 

35 

03 

CO 

35 

6Z 

26 

149 

26 

70 

00 

00 

26 

79 

6Z 

70 

CM 

1 ^ 

CM 

CD 

46 

143 


- 



03 

CO 


CO 

00 

o 

CO 

03 

03 

00 

CO 

Ojunp 

53 

149 

491 

''t 

03 

00 

26 

53 

26 

96 

26 

70 

79 

26 

CO 

70 

26 

LO 

1^ 

CD 


56 

LO 

- 

CM 

00 

CD 

CO 

- 

CO 

00 

03 

CO 


00 

00 

O 

CM 

oAbiai 

79 

70 

272 

35 

00 

35 

53 

44 

114 

53 

44 

140 

140 

35 

79 

350 

35 

00 

00 

03 

00 

CO 


CM 


CD 

LO 

CO 

CD 

LO 

CD 

CD 


03 

40 


O 

IMQV 

140 

526 

902 

53 

00 

00 

70 

00 

00 

00 

00 

00 

00 

149 

96 

CO 

CD 

CM 

CO 

CD 

CM 

70 

00 

00 


79 

228 

CD 

60 

103 

CD 

O 

00 

O 

o 

o 


- 

30 

30 

00 

O 


03 

26 

OZJBIAI 

114 

289 

622 

105 

96 

105 

96 

CD 

CD 

96 

140 

00 

00 

324 

412 

70 

105 


70 

254 

00 

33 


CM 

- 

CM 

- 

03 

- 

CD 

o 

37 

47 

00 

CM 


00 

29 

oj0jq0d 

149 

876 

2111 

70 

140 

403 

123 

263 

114 

140 

6Z 

350 

254 

70 

210 


105 

894 


100 

241 

00 

CD 

46 


30 

CO 

CD 

03 

40 

29 

00 

24 


CM 

ZOl 

OJ 0 U 3 

149 

447 

1428 

53 

123 

456 

105 

307 

79 

175 

6Z 

158 

350 

53 

105 


114 

569 


5 

163 

CD 

^f- 

52 

CM 

35 

03 

20 

03 

00 

40 

CD 

CM 


00 

65 

Estacion 

CM 

03 

CO 

CD 

LO 

00 

CD 

O 

_Q 

00 

03 

00 

00 

00 


CO 

CM 

E 

CD 

E 

LO 

c 

o 
c n 

CM 

C/3 

D 

00 

> 

CO 

> 

xt" 

N 
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Tabla 5-15 Velocidad promedio del viento en m/s obtenidas en el proyecto SWERA. 


Estacion 

Enero 

Febrero 

Marzo 

Abril 

Mayo 

Junio 

Julio 

Agosto 

Septiembre 

Octubre 

Noviembre 

Diciembre 

Promedio Anual 

a12 

2.39 

2.27 

2.19 

2.37 

1.83 

1.62 

1.73 

1.72 

1.58 

1.82 

1.99 

2.41 

1.95 

a31 

2.8 

3.83 

2.4 

2.83 

1.64 

1.71 

2.92 

2.65 

1.77 

3 

3.57 

3.61 

2.65 

a35 

4.55 

4.76 

3.09 

3.55 

2.35 

2.67 

5 

3.95 

3.55 

4.28 

5.08 

5.38 

3.91 

blO 

1.65 

1.76 

1.78 

1.6 

1.37 

1.24 

1.41 

1.41 

1.31 


1.42 

1.7 

1.12 

g3 

1.47 

1.55 

1.45 

1.45 

0.95 

0.78 

0.75 

0.94 

0.8 

0.86 

0.97 

1.42 

1.12 

hi 3 

2.96 

1.85 

1.76 

1.65 

1.36 

1.26 

1.72 

1.76 

1.49 

2.17 

2.43 

1.99 

1.74 

h8 

1.87 

1.86 

1.81 

1.74 

1.43 

1.19 

1.36 

1.22 

1.11 

1.28 

1.66 

1.79 

1.55 

118 

2.76 

2.62 

1.87 

1.48 

1.25 

1.45 

1.85 

1.82 

1.65 

2.01 

2.96 

3.11 

2.08 

14 

1.51 

1.76 

1.69 

1.63 

1.68 

1.06 

1.09 

1.09 

1.07 

0.09 

1.24 

1.48 

1.34 

m23 

2.16 

2 

2 

1.97 

1.69 

1.9 

2.16 

2.13 

1.83 

2.12 

2.12 

2.79 

2.12 

m6 

1.6 

1.55 

1.69 

1.64 

1.26 

0.98 

1.09 

1.1 

0.87 

0.83 

0.89 

1.39 

1.26 

n15 

2.38 

3.34 

3.12 

2.67 

2 

1.7 

1.65 

1.51 

1.43 

1.39 

1.57 

2.16 

1.87 

slO 

3.12 

2.84 

3.06 

2.81 

2.13 

1.89 

2.1 

2.05 

1.87 

2.11 

2.67 

3.11 

2.44 

s27 

1.66 

1.73 

1.68 

1.71 

1.24 

1.28 

1.28 

1.29 

1.23 

1.31 

1.49 

1.64 

1.5 

ull 

1.79 

1.87 

1.74 

1.77 

1.67 

1.62 

1.73 

1.52 

1.48 

1.49 

1.57 

1.65 

1.64 

v13 





2.41 

1.44 

1.35 

1.53 

1.29 

1.17 

— 

— 

1.5 

v6 

1.86 

1.8 

1.49 

1.64 

1.18 

1.18 

1.62 

1.26 

1.13 

1.06 

1.34 

1.72 

1.43 

z4 

3.54 

4.54 

1.5 

2.66 

1.85 

2.47 

3.14 

3.21 

2.26 

3.28 

3.79 

4.34 

2.88 


Tabla 5-16 Datos de viento (W/m 2 y KWh/m 2 ) para cada una de las estaciones climaticas del pais. 
Calculos basados en la tabla 1-1 (Duracion aproximada de las estaciones climaticas. 


Estacion 

meteorologica 

Seca 

Seca - 
Lluviosa 

Lluviosa 

Lluviosa 

Seca 

a12 

15 

131 

10 

88 

6 

53 

10 

88 

a31 

66 

578 

11 

96 

25 

219 

74 

648 

a35 

170 

1489 

29 

254 

88 

771 

186 

1629 

blO 

7 

61 

5 

44 

6 

53 

5 

44 

93 

13 

114 

3 

26 

5 

44 

10 

88 

hi 3 

32 

280 

4 

35 

7 

61 

39 

342 

h8 

12 

105 

6 

53 

3 

26 

7 

61 

118 

31 

272 

5 

44 

9 

79 

40 

350 

14 

8 

70 

14 

123 

3 

26 

5 

44 

m23 

23 

201 

5 

44 

13 

114 

— 

— 

m6 

9 

79 

6 

53 

3 

26 

2 

18 

n15 

24 

210 

23 

201 

9 

79 

4 

35 

slO 

37 

324 

18 

158 

10 

88 

19 

166 

s27 

7 

61 

7 

61 

1 

9 

4 

35 
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Estacion 

meteorologica 

Seca 

Seca - 
Lluviosa 

Lluviosa 

Lluviosa 

Seca 

ull 

12 

105 

10 

88 

6 

53 

5 

44 

v13 

— 

— 

54 

473 

12 

105 

4 

35 

v6 

10 

88 

5 

44 

5 

44 

4 

35 

z4 

52 

456 

13 

114 

22 

193 

65 

569 



fit 


swera 


Tabla 5-17 Datos de velocidad del viento en m/s para cada una de las estaciones climaticas del pais. 

Calculos basados en la tabla 1-1. 


Estacion 

meteorologica 

Seca 

Seca - 
Lluviosa 

Lluviosa 

Lluviosa 

Seca 

a12 

2.29 

1.94 

1.68 

2 

a31 

3.12 

1.78 

2.17 

3.93 

a35 

4.31 

2.41 

3.8 

4.57 

blO 

1.67 

1.58 

1.45 

1.45 

93 

1.44 

1.11 

0.91 

1.4 

hi 3 

2.23 

1.4 

1.6 

2.67 

h8 

1.81 

1.46 

1.25 

1.63 

118 

2.5 

1.19 

1.7 

2.71 

14 

1.34 

1.78 

1.07 

1.29 

m23 

2.38 

1.64 

1.96 

— 

m6 

1.52 

1.32 

1 

0.86 

n15 

2.36 

2.28 

2.37 

1.37 

slO 

3.02 

2.18 

1.97 

2.42 

s27 

1.66 

1.53 

0.73 

1.47 

ull 

1.72 

1.68 

1.58 

1.5 

v13 

— 

2.81 

1.44 

1.2 

v6 

1.67 

1.26 

1.37 

1.27 

z4 

3.35 

2.06 

2.66 

3.79 
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Tabla 5-18 Tabla comparativa de datos de viento m/s, W/m 2 , KWh/m 2 para diferentes proyectos. 


Estacion 

Promedio anual SWERA 

Promedio anual UPES 

Promedio anual MINEC 3 

a12 

1.95 

12 

105 

2.13 

31 

280 

2.1 

31.8 

278 

a31 

2.65 

45 

394 

3.14 

70 

629 

3.14 

65.4 

573 

a35 

3.91 

123 

1077 

4.38 

116 

1040 

4.33 

118.7 

1040 

blO 

1.12 

8 

70 

1.4 

15 

133 

1.41 

15.2 

133 

93 

1.12 

8 

70 

1.41 

16 

147 

1.41 

16.6 

146 

hi 3 

1.74 

19 

166 

2.06 

31 

275 

2.03 

31.4 

275 

h8 

1.55 

8 

70 

2.14 

31 

277 

2.16 

32 

280 

118 

2.08 

22 

193 

2.24 

36 

319 

2.26 

36.7 

319 

14 

1.34 

6 

53 

1.5 

17 

155 

1.5 

17.6 

154 

m23 

2.12 

17 

149 

2.42 

37 

334 

2.42 

39.7 

348 

m6 

1.26 

6 

53 

1.5 

17 

152 

1.54 

17.5 

153 

n15 

1.87 

16 

140 

2.55 

45 

407 

2.55 

46.2 

405 

slO 

2.44 

22 

193 

2.16 

32 

289 

2.17 

33 

289 

s27 

1.5 

5 

44 

1.55 

18 

159 

— 

— 

— 

s5 

— 

— 

— 

— 

— 

— 

1.98 

27.7 

243 

ull 

1.64 

8 

70 

1.94 

26 

230 

1.95 

26.2 

230 

v13 

1.5 

14 

123 

2.62 

48 

428 

2.62 

48.7 

426 

v6 

1.43 

7 

61 

1.49 

17 

149 

1.49 

17.1 

150 

z4 

2.88 

33 

289 

2.87 

57 

510 

2.81 

57.9 

508 

t6 

— 

— 

— 

2.87 

53 

475 

2.87 

55.9 

482 

pi 1 

— 

— 

— 

2.23 

34 

300 

2.23 

34.3 

300 


Tabla 5-19 Serie de datos para mediciones eolicas. Los datos en azul son los datos mensuales de la 
base de datos UPES y los datos en rosado son los datos horarios recuperados en el proyecto SWERA- 

UCA. 


Santa Ana El Palmar 




o 

■ Tornados del CD elaborado por el consultor Humberto Rodriguez para el proyecto del MINEC “Implementation de 
medidas para el desarrollo de un mercado national sostenible para servicios de energia en el area rural”. 2005. 
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Planes de Montecristo 



Anos de registros 



Meses 


Cerro Verde 



Anos de registros 



Meses 


Cerron Grande 



Anos de registros 



Nueva Concepcion 



Anos de registros 



Apaneca Santa Leticia 



Anos de registros 



Meses 
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Ahuachapan 



Anos de registros 



Meses 


El Boqueron 



Anos de registros 



2.0 
1.5 
1.0 
0.5 

0 2 4 6 8 10 12 

Meses 



San Andres 



Anos de registros 



2.0 
1.5 
1.0 
0.5 
0.0 

0 2 4 6 8 10 12 

Meses 



Cerro Cacahuatique 



Anos de registros 



2.5 
2.3 
2.1 
1.9 
1.7 

1.5 

0 2 4 6 8 10 12 

Meses 



Beneficio el Papalon 




Anos de registros Meses 
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La Union 


I 



I 



Anos de registros 


Meses 


Aeropuerto llopango 


I 



i 


Anos de registros 



Estacion Matriz 


I 



I 



Anos de registros 


Meses 


Beneficio La Carrera 


I 



I 


Anos de registros 



Meses 


Apastepeque 


I 



I 


Anos de registros 



Meses 
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Santa Cruz Portillo 


I 



i 


Anosde registros 



La Galera 


I 



I 


A nos de registros 



Meses 
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Mapa Eolico de El Salvador (NREL 2004) 

Potencia media anual a 50 m.a.g.l. 
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Mapa Eolico de El Salvador (NREL 2004) 

Potencia media anual a 50 m.a.g.l. 
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6. Conclusiones y recomendaciones. 

1 . Para la estimacion de los mapas de radiacion solar se mejoro la serie de datos con mediciones 
de 5 a 14 anos del proyecto de Solarimetria a una nueva serie de datos de 5 a 30 anos con el 
proyecto SWERA. Para los datos de luz solar se utilizo la serie del proyecto de solarimetria de 5 
a 25 anos completado con dos anos de datos (para cada una de las estaciones) digitalizados con 
el proyecto SWERA. 

2. La tecnologia de SIG permitio que se hicieran de mejor manera las asignaciones de constantes 
de Angstrom a estaciones con mediciones de luz solar unicamente. Para ello se utilizaron 
diferentes criterios como: el mapa de zonas de vida de Holdrige, las elevaciones de las 
estaciones (altitud) y las cercanias entre las estaciones. Para la generacion de los mapas de 
potencial solar (luz solar y radiacion solar) se utilizo la informacion de 22 estaciones 
meteorologicas por medio de la interpolacion de Kriging. 

3. Se ha realizado un esfuerzo importante en la recuperacion de los datos digitales de sol y viento. 
El inventario de radiacion solar esta completo y falta terminar de digitalizar las series de datos de 
luz solar y bandas de viento. En el documento se detallan los porcentajes de avances en el 
proceso de digitalizacion. 

4. Para la determinacion del potencial eolico del pais se utilizo el software WAsP , el cual permitio 
modelar la climatologia regional alrededor de tres estaciones meteorologicas identificadas como 
de mayor potencial eolico: cerro verde, planes de Montecristo y la galera. Ademas, se uso esta 
herramienta de software para calcular los atlas de viento de las 18 estaciones meteorologicas 
para diferentes meses del ano y estaciones climaticas. 

5. Se desarrollo un sistema informatico que permite incorporar datos nuevos de sol y viento a una 
base de datos la cual puede ser consultada, generar diferentes tipos de reportes, etc todo esto 
en un ambiente de red lo que permite aumentar la productividad de los encargados de mantener 
actualizadas las bases de datos solar y eolica del pais. Tambien este sistema podria ser utilizado 
por otros paises de Centroamerica que utilizan instrumentos de medicion similar a los de nuestro 
pais. 

6. Los mapas generados de brillo solar y radiacion solar deben analizarse dentro del contexto que 
representan modelos de superficies generados por el metodo de interpolacion de Kriging. Otros 
metodos de interpolacion de seguro produciran modelos de superficies diferentes a los 
generados en este estudio. Tambien debe considerarse que estos metodos geoestadisticos 
dependen del criterio del investigador para ajustar el semivariograma experimental a los modelos 
tradicionales de semivariogramas (circular, estable, esferico, exponencial, gaussiano, etc.). 

7. Las bandas escaneadas para el caso de la radiacion y el viento constituyen una importante base 
de datos en digital que permitiria en un futuro eliminar fisicamente los inventarios en papel de 
estas imagenes. 

8. Para cada uno de los parametros de estudio de sol y viento se establece un procedimiento 
informatico para su recuperacion lo cual permite que la division de meteorologia del SNET 
actualice continuamente estos parametros. 

9. La suma de las constantes de Angstrom (a y b) encontradas en este estudio andan alrededor de 
0.65, lo que esta de acuerdo a lo que establece la literatura para climas tropicales y humedos 
como el de nuestro pais. 

10. El analisis de la climatologia regional nos permite concluir que el mapa eolico elaborado por 
NREL en el ano 2004 para el pais es un excelente instrumento para la determinacion de sitios 
adecuados para la explotacion eolico tanto a nivel residencial como de granjas eolicas. 

1 1 . En los tres lugares en que se determino la climatologia regional se encontraron sitios con valores 
altos de potencia y velocidad que no se evidenciaban con los valores registrados por las 
estaciones meteorologicas por lo que este estudio no descarta la posibilidad de la creacion de 
granjas eolicas en El Salvador. 

12. El mapa eolico de NREL es una excelente gufa para la ubicacion de estaciones meteorologicas 
con el objetivo de afinar las mediciones de viento en las zonas identificadas en el mapa eolico. 

13. Se recomienda que la divisiones de meteorologia del SNET recupere la informacion horaria de 
viento utilizados por los estudios de UPES (1987) y MINEC (2005) de tal forma de crear una sola 


Proyecto SWERA-EL SALVADOR 


6-1 



14. 


15. 


16. 


17 . 


w 


Hnb»Anfrn2r:i 

HrvomNilukn 


swera 


base de datos digital en viento que permita hacer mejores estimaciones que las realizadas oir el 
proyecto SWERA. 

Comparando las series de datos de viento del proyecto SWERA con las historicas mostradas en 
otros estudios, indican un deterioro de los instrumentos de medicion (anemdmetros) en algunas 
de las estaciones, lo cual hace necesario que se tome en cuenta al utilizarlos los datos en futuras 
investigaciones. 

El modelaje eolico desarrollado en este estudio no considero perdidas por obstaculos alrededor 
de las estaciones meteorologicas y estas debieron haber sido consideradas al observar las rosas 
de vientos de algunas estaciones que lo evidenciaron. Esto dificulto hacer un detalle completo de 
las huellas meteorologicas de las estaciones del pais con registros de viento. 

La incorporacion de lenguaje informatico SQL de bases de datos relacionales y herramientas de 
sistemas de informacion geograficas simplified y sistematizo el manejo y calculos de los datos 
meteorologicos presentados en este estudio. Se sistematizo adicionalmente los formatos de 
entrada y salida de WAsP para su representacion y analisis con sistemas de informacion 
geografica. 

El grupo de investigadores de la UCA agradece a la division de meteorologia de SNET (MARN) 
por el acceso a los inventarios de informacion solar y eolica de el pais y el apoyo tecnico 
oportuno. Asf mismo se agradece al proyecto SWERA por el fortalecimiento de las capacidades y 
habilidades adquiridas por los investigadores en los temas solar y eolico. 


Proyecto SWERA-EL SALVADOR 


6-2 




7. Bibliografia. 


1. Burrough, P.A., “Principles of Geographic Information Systems”, Oxford University Press, 
1998. 

2. Manuales de ArcView y Spatial Analyst, ESRI. 

3. Manuales de Idrisi Kilimanjaro, 2003. 

4. PanatierY., VARIOWIN: Software for Spatial Data Analysis in 2D, Springer, 1996. 

5. Informe del proyecto de Solarimetria publicado por la UES en septiembre de 1989. 
Financiamiento de OEA. Ejecucion del Departamento de Fisica de la Universidad de El 
Salvador. 

6. Duffie, John, Beckman, W. “Solar engineering of thermal processes”. John Wiley & Sons. 
1980. 

7. Modulos de radiacion solar publicados en el sitio web 
http://www.courses.ait.ac.th/ED06.22/course1/lecs/ 

8. lb Troen y Erick Lundtang Petersen, “El atlas eolico europeo”, Meteorology and Wind Energy 
Department, Riscp National Laboratory, 1990. 

9. WAsP courses and manuals, “Wind Atlas Analysis and Application Program”, course notes 
2002 y 2004. 

10. “Small wind electric systems: an Ohio consumer's guide”, US Department of Energy, 2005. 

11. Eduardo Morales Manchame, “El viento en Centroamerica”, principios de los anos 80. 

12. Mick Sagrillo, “Site Analisis for Wind Generators”, abril-mayo de 1994. 

13. Alexandre de Lemos Pereira, “Site selection for wind power project in El Salvador”, Noviembre 
de 1994, convenio UCA-CEL-MARN. 

14. Humberto Rodriguez, “Implementacion de medidas para el desarrollo de un Mercado nacional 
sostenible para servicios de energia en el area rural de El Salvador”, Ministerio de Econorma, 
Mayo de 2005, CD de divulgacion. 

15. Seccion Diagnostico-Clima, vientos en El Salvador, “Plan Nacional de ordenamiento y 
desarrollo territorial”, 2005. 

16. Lisa M. Daniela, Sara Johnson, “Harvest the wind: A wind energy handbook for Illinois”, por 
WINDUSTRY, 2004. 

17. Trabajo de graduacion “Extraccion de agua subterranea con energia eolica en El Salvador” 
realizado por estudiantes de la Universidad Politecnica de El Salvador (UPES) en 1987. 

18. “Primer plan nacional de desarrollo energetico integrado 1988-2000”, publicado en enero de 
1988 por la Comision Ejecutiva Hidroelectrica del Rio Lempa (CEL). 


Proyecto SWERA-EL SALVADOR 


7-1 



